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高 速 飞行 器 空气 动力 学 数值 分 析 
摘要 

高 速 飞行 器 空气 动力 学 数值 分 析 是 一 种 以 计算 流体 力学 为 基础 的 准确 、 高效 和 实用 
的 方法 。 目 前 国内 外 对 计算 流体 力学 的 研究 有 从 求解 欧 拉 方 程 到 求解 NS 方程 ， 从 求解 
二 维 问题 到 求解 三 维 问题 的 趋势 ， 它 要 求 数值 求解 NS 方程 组 来 模拟 粘性 效应 。 

本 文 基于 高 速 飞 行 器 空气 动力 学 数值 分 析 ， 选 用 了 雷诺 平均 NS 方程 (RANS)， 采 用 
k-c 模型 ， 定 义 专门 的 边界 模型 条 件 ， 以 压力 、 速 度 、 温 度 、k 和 为 基本 量 进 行 迭代 
求解 ， 对 松弛 因子 的 取 值 范围 合理 缩小 化 ， 根 据 SIMPLE 算法 和 求解 模型 研制 了 相应 的 
软件 。 对 高 速 飞行 器 的 流 场 进行 数值 分 析 ， 得 出 了 在 0 攻 角 下 ， 两 种 高 速 飞行 器 流 场 的 
速度 分 布 、 压 力 分 布 以 及 温度 分 布 ， 并 对 来 流 为 不 同 马 赫 数 时 ， 两 种 高 速 飞 行 器 数值 分 
析 的 阻力 系数 和 风 洞 实验 阻力 系数 进行 比较 ， 比 较 结果 表明 有 一 致 性 ， 从 而 验证 了 模型 
及 软件 的 正确 性 。 在 此 过 程 中 主要 做 的 以 下 几 方面 的 研究 : 

L 给 出 流体 数学 模型 ， 定 义 专 门 的 初始 边界 模型 条 件 ; 

2， 各 种 满 流 模型 的 对 比 ; 

3， 详 细 介 绍 了 三 维 SIMPLE 算法 , 给 出 求解 器 的 设计 方法 ; 

4. 求解 了 两 种 高 速 飞行 器 气动 参数 ， 和 风 洞 实验 的 气动 参数 作 了 比较 ， 结 果 一 致 
性 比较 好 。 

关键 词 : 计算 流体 力学 , 高 速 飞 行 器 ，SIMPLE 


High Speed Aircraft Aerodynamics Numeric Analysis 
ABSTRCT 


High speed aircraft aerodynamics numeric analysis is a kind of exact effective and 
applied method which based on CFD. Nowadays ,The trend of the CFD research domain is 
from solving Euler equation to solving NS equation, from 2 dimensions to 3 dimensions, 
which requires numerically solving NS equation group to simulate viscosity domino effect. 

Based on High Speed Aircraft Aerodynamics Numeric Analysis , RANS equation is 
chose, *-€ equation model is adopted, special model boundary conditions are defined, 
pressure, velocity, temperature, k and £ as basic quantity are used to solve it by iteration, 
relaxation factor is downsized legally, corresponding software is designed to simulate the air 
field of the high speed aircraft according to SIMPLE algorithm and solving model, velocity 
field ,pressure field and the temperature field in the flow field around high speed aircraft is 
gained when attack angle is 0 degree and comparison is made on drag coefficient between 
numerical analysis and wind tunnel test when attack angle is 0 degree and flow comes in 
different velocity. Finally the consistency is got and the legitimacy of the model and software 
are checked .I did the research on the following jobs: 

1. giving mathematic model, defining special model boundary conditions; 

2. Comparing with all kinds of turbulence models and giving the design method of 
solver; 

3. Introducing SIMPLE algorithm detailedly; 

4, Solving pneumatic parameter of high speed aircraft and making comparison on drag 
coefficient between numerical analysis and wind tunnel test, finally getting the 
consistency. 


KEYWORD: CED, high speed aircraft, SIMPLE 
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第 一 章 绪论 


高 速 飞行 器 空气 动力 学 数值 分 析 是 一 种 以 计算 流体 力学 为 基础 的 准确 、 高 效 和 实用 
的 方法 。 在 工程 应 用 的 需求 和 推动 下 ， 计 算 流 体力 学 随 着 计算 机 的 飞速 发 展 和 计算 格式 
的 不 断 更 新 而 得 到 迅 狐 发 展 ， 同 时 在 工程 实际 中 的 应 用 范围 也 越 来 越 广泛 ， 逐 步 与 地 
面试 验 、 飞 行 实验 一 起 成 为 飞行 器 性 能 分 析 、 飞 行 器 空气 动力 学 设计 和 飞行 器 设计 的 三 
大 工具 。 应 用 计算 流体 力学 对 高 速 飞行 器 空气 动力 学 数值 分 析 有 重要 意义 ， 

首先 是 经 济 、 快 捷 。 计 算 流体 力学 (CPD) 应 用 到 航空 领域 具体 为 计算 空气 动力 学 ， 
它 既 有 助 于 我 们 理解 复杂 流动 的 物理 本 质 ， 同 时 也 是 对 风 洞 实验 的 一 个 补充 ， 可 以 大 大 
缩短 飞行 器 设计 和 研制 的 周期 ， 降 低 整 个 飞行 器 的 设计 费用 。 

其 次 是 可 直接 获取 基本 方程 的 解 。 能 给 出 流 场 包括 对 绕 流 物体 的 力 、 热 、 密 度 等 所 
有 详尽 的 数据 , 能 够 方便 地 识别 一 些 关 键 参 数 的 影响 和 探索 力学 现象 相互 作用 的 结果 和 和 
规律 。 在 概念 设计 、 初步 设计 阶段 , 用 计算 流体 动力 学 可 以 较 快 地 进行 技术 可 行 性 分 析 ， 
多 种 方案 的 筛选 。 这 一 阶段 如 果 进 行 试 验 ， 在 时 间 和 经 费 方向 都 将 付出 很 大 代价 ， 而 且 
不 可 能 像 计 算 那 样 在 广阔 的 范围 内 筛选 ， 通 过 计算 ， 在 广阔 的 范围 内 筛选 后 再 用 试验 进 
行 精 选 ， 可 以 更 有 效 她 保证 质量 和 和 避免 技术 风险 。 

总 之 ， 在 高 速 飞行 器 动力 学 数值 分 析 中 ,广泛 使 用 计算 流体 动力 学 不 仅 可 以 保证 质 
量 ， 减 少 和 避免 风险 、 而 且 可 能 以 较 小 的 代价 获取 飞行 器 的 优化 。 


1.2 国内 外 研究 现状 
1. 2. 1 计算 流体 力学 的 研究 
随 着 现代 航空 航天 技术 的 高 速 发 展 , 实际 飞行 中 的 飞行 器 周围 的 流 场 变 得 越 来 越 复 


杂 。 复 杂 流 动 现象 的 主要 研究 手段 是 实验 分 析 法 和 数值 分 析 。 实 验 分 析 法 的 结果 比较 真 
实 ， 是 检验 理论 结果 和 计算 结果 的 重要 标准 ， 但 存在 着 耗资 昂贵 、 周 期 长 、 仿 真 条 件 受 
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限 、 不 能 显示 和 测量 所 有 的 流动 现象 、 不 能 同时 满足 几 个 相似 准则 、 有 测量 误差 等 不 足 
之 处 。 而 对 于 描述 复杂 流动 现象 的 数学 模型 方程 ， 用 解析 方法 求解 是 不 可 能 的 。 因 此 ， 
在 六 十 年 代 产 生 了 用 数值 方法 模拟 和 分 析 流 动 现 象 的 学 科 : 计算 流体 力学 (CFD) 的 发 展 
一 般 可 以 分 为 4 个 阶段 ” COO 线性 化 无 粘 近似 ; (2) 非 线性 无 粘 近 似 ， 如 欧 拉 方程 ; 
(3) 雷诺 平均 RANS 方 程 ，《〈4) 非 定常 N-S 方 程 。 例 如 六 十 年 代 研究 了 线性 化 方程 的 求 
解 ， 七 十 年 代 进 行 了 无 粘 位 势 流 动 方程 的 求解 ， 八 十 年 代 开 始 对 Euler/N-S 方 程 的 求解 ， 
同时 数值 模拟 方法 在 工程 实践 中 得 到 广泛 发 展 。 

近 几 年 来 ，CFD 技 术 发 展 主要 表现 在 以 下 几 个 方面 2 s 

CL) 从 求解 欧 拉 方程 发 展 到 Navier-Stokes 方 程 ;从 求解 二 维 问题 发 展 到 三 维 问题 ; 
从 简单 翼 型 、 机 辟 发 展 到 全 机 实用 外 形 ， 并 且 能 模拟 军用 飞机 与 外 挂 导弹 及 民用 飞机 与 
短 舱 之 间 相 互 干扰 的 复杂 流 场 ; 从 计算 小 迎 角 发 展 到 大 迎 角 ， 从 定常 流 场 数值 模拟 发 展 
到 复杂 的 非 定常 流 场 数值 模拟 和 相关 问题 求解 ;从 亚 、 超 音速 发 展 到 处 理 跨 音速 、 高 超 
音速 问题 。 

(2) 从 空间 、 时 间 的 一 阶 精度 计算 发 展 更 高 阶 精度 计算 ， 

(3) 大 部 分 分 析 流 体 软件 都 采用 有 限 体积 的 离散 方法 ， 目 前 流行 的 大 型 商业 化 CFD 
软件 中 只 有 FIDAP 是 使 用 有 限 元 法 的 ， 而 CFX、PHOENICS、STAR 一 CD、FLUENT 等 均 采 用 有 
限 容积 法 

(4 从 利用 结构 网 格 发 展 到 利用 非 结构 网 格 进行 流 场 离散 ;网 格 形态 从 规则 的 四 边 
形 、 六 面体 结构 网 格 已 发 展 到 更 具有 灵活 性 的 三 角形 、 由 面体 、 多 边 形 、 多 面体 等 非 结 
构 化 网 格 ， 更 具有 对 复杂 接口 的 普遍 适用 性 。 

(5) 网 格 可 以 分 成 许多 不 同 结构 的 块 (B10CK) ,允许 不 同 结构 块 之 则 的 数据 通讯 及 交 
换 ， 从 而 可 以 将 复杂 外 形 进行 规则 的 划分 ， 以 适用 于 复杂 的 几何 形状 。 

(6) 网 格 单元 间距 的 划分 已 实现 了 自 适应 化 ， 可 以 根据 问题 物理 量 的 变化 梯度 优化 
分 布 网 格 ， 提 高 解 题 精度 及 效率 。 并 可 以 进行 局 部 加 密 细 化 ， 提 高 局 部 的 分 辩 率 ， 更 好 
地 捕捉 流动 现象 及 机 理 。 

(7) 在 间断 面 、 线 附近 或 存在 计算 不 稳定 情形 下 使 用 自动 人 工 粘 性 技术 ， 有 利于 获 
得 平滑 、 稳 定 的 数值 解 。 

(8) 采 用 多 重 网 格 加 速 来 提高 收敛 的 速度 及 稳定 性 。 
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(9) 求 解 的 技术 也 从 单机 运算 发 展 到 多 CPU 工作 站 ， 分 布 式微 机 网 络 ， 大 规模 机 群 并 
行 计 算 ， 新 的 CFD 程 序 设计 思想 逐渐 实用 化 。 其 中 包括 面向 对 象 的 设计 思想 和 各 种 并 行 
计算 设计 方法 ， 如 : 分 区 并 行 计算 、 共 享 内 存 的 多 CPU 并 行 计算 、 采 用 分 布 式 内存 的 大 
规模 并 行 计 算 、 采 用 工作 站 集群 方式 进行 的 网 络 化 大 规模 并 行 计 算 等 ， 代 表 了 当今 发 展 
的 潮流 。 


1.2.2 ” 弹 箭 空气 动力 学 的 研究 现状 


自 二 十 世纪 六 十 年 代 以 来 CFD 技 术 得 到 飞速 发 展 , 其 原动力 是 不 断 增 长 的 工业 需求 ， 
而 航空 航天 工业 自始至终 是 最 强大 的 推动 力 。 传 统 飞行 器 设计 方法 试验 昂贵 、 费 时 ， 所 
获 信息 有 限 ， 追 使 人 们 需要 用 先进 的 计算 机 模拟 手段 指导 设计 ， 大 量 减 少 原型 机 试验 ， 
缩短 研发 周期 ， 节 约 研 究 经 费 。 四 十 年 来 ，CFD 在 满 流 模 型 、 网 格 技术 、 数 值 算 法 、 可 
视 化 、 并 行 计算 等 方面 取得 飞速 发 展 ， 并 给 工业 界 带 来 了 革命 性 的 变化 。 如 国外 飞机 厂 
商用 CFD 取 代 大 量 实 物 试验 ,如 美国 战斗 机 YF-23 采 用 CFD 进 行 气动 设计 后 比 前 一 代 YF-17 
减少 了 60% 的 风 涧 试验 量 。 目 前 在 航空 、 航 天 、 汽 车 等 工业 领域 ， 利 用 CFD 进 行 的 反复 设 
计 、 分 析 、 优 化 已 成 为 标准 的 必 经 步骤 和 手段 “*” 。 

长 期 以 来 ， 飞 行 器 的 气动 弹性 问题 都 是 飞行 器 设计 必须 考虑 的 重要 课题 。 从 固定 机 
要 飞行 器 问世 的 第 一 天 起 ， 人 类 就 开始 面临 飞 行 器 的 气动 弹性 问题 , 1903 年 Langley 的 单 
田 机 在 作 有 动力 的 首次 飞行 时 , 就 因为 典型 的 静 气 动弹 性 扭转 发 散 而 折断 了 机 姻 。 近 100 
年 来 ， 随 着 飞行 器 研制 的 不 断 进 步 ， 飞 机 的 气动 弹性 问题 也 日 益 成 为 影响 飞行 器 飞行 性 
能 的 难题 ， 例 如 美国 的 F18-EIF 战 斗 机 的 研制 过 程 中 ， 机 杜 在 一 定 迎 角 下 飞行 时 就 遇 到 
静 气 动弹 性 扭转 发 散 的 问题 。 我国 在 国产 飞机 的 自行 研制 和 型 号 改进 的 时 候 ， 也 不 只 一 
次 的 外 到 过 拌 振 ， 喻 鸣 等 气动 弹性 问题 。 由 于 飞行 器 的 气动 弹性 问题 常常 严重 影响 飞行 
器 的 飞行 性 能 和 乘坐 舒适 性 ， 严 重 的 还 会 造成 飞机 结构 破坏 ， 威 胁 乘 员 生命 安全 ， 所 以 
必须 深入 研究 ， 认 真 解决 。 目 前 ,我 国正 在 预 研 的 第 四 代 战 斗 机 要 求 具 有 跨 音 速 范围 的 
高 机 动 性 和 超 音 速 巡 航 特性 , 正在 研制 的 大 型 运输 机 、 训 炸 机 、 民用 客机 和 新 型 无 人 机 ， 
一 般 都 具有 大 展 弦 比 或 超大 展 弦 比 机 由 ， 使 得 机 村 的 结构 弹性 变形 较 大 , 大 量 轻 质 复合 
材料 的 广泛 应 用 ， 使 得 飞行 器 的 结构 刚度 特性 发 生变 化 。 这 些 新 情况 的 出 现 使 飞行 器 设 
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计 中 气动 弹性 问题 日 益 突出 。 以 往 以 刚性 飞机 为 设计 对 象 加 以 弹性 修正 的 设计 方法 不 能 
满足 新 型 飞机 的 设计 要 求 ， 需 要 以 真实 的 弹性 飞行 器 为 对 象 直接 进行 飞行 器 气动 特性 研 
究 分 析 ， 这 些 急 待 解决 的 新 课题 就 使 本 文 的 研究 具有 了 非常 紧迫 的 现实 意义 。 

由 于 飞行 器 气动 弹性 力学 属于 多 学 科 交 叉 的 学 科 , 其 研究 的 问题 同时 涉及 到 流体 力 
学 、 国 体力 学 以 及 控制 理论 等 多 个 学 科 的 难点 ， 因 此 飞行 器 的 气动 弹性 一 直 以 来 都 是 飞 
行 器 设计 工作 者 必须 面 对 的 重要 研究 难题 。 在 以 往 的 飞行 器 设计 过 程 中 ， 飞 行 器 设计 工 
作者 主要 是 通过 实验 来 获取 飞行 器 气动 弹性 特性 。 由 于 气动 弹性 的 风 洞 实验 较 之 定常 流 
场 实验 ， 需 要 设计 专门 的 实验 仪器 ， 实 验 过 程 非常 复杂 ， 精 度 也 有 限 ， 而 飞行 器 飞行 试 
验 危险 系数 很 大 。 因 此 飞行 器 设计 工作 者 一 般 很 难 准确 得 到 飞行 器 的 气动 弹性 特性 , 只 
能 主要 通过 一 些 工 程 计 算 方法 ， 从 结构 力学 角度 对 飞行 器 的 气动 弹性 特性 进行 预测 。 例 
如 目前 大 多 数 工程 上 应 用 的 频 域 方法 ， 其 气动 力 的 求解 都 是 基于 线性 的 升力 面 理论 ， 气 
动力 的 求解 只 能 满足 亚 音速 ， 低 马赫 数 的 气动 弹性 计算 。 而 在 飞行 器 跨 音 速 阶段 ， 由 于 
率 流 附 面 层 ， 跨 音速 流动 的 分 离 ， 激 波 的 不 稳定 波动 ， 以 及 激 波 附 面 层 的 于 扰 等 现象 的 
出 现 ， 气 动力 将 会 出 现 复杂 的 非 线性 特性 ,气动 力 的 非 线性 变化 使 得 传统 的 工程 算法 求 
解 的 气动 力 与 实际 情况 相差 甚 远 。 例 如 额 振 , RS TOME STEPS EFE BLUR 
现 特有 的 “四 坑 ” 现 象 ， 使 飞行 器 在 跨 音 速 阶段 发 生 颤 振 的 可 能 性 大 大 增加 。 过 去 飞行 
器 设计 人 员 在 设计 阶段 由 于 很 难 准确 获得 飞行 器 跨 音 速 的 气动 弹性 特 人 性， 只 能 通过 飞行 
器 亚 音速 气动 弹性 特性 进行 估 测 ， 严 格 控制 飞行 包 线 来 避免 跨 音 速 的 气动 弹性 不 稳定 。 
当 飞 行 器 上 天 后 遇 到 气动 弹性 问题 ， 也 只 能 通过 一 些 经 验 办 法 ， 如 增加 机 辟 刚 度 ， 加 扰 
流 片 等 来 加 以 控制 。 这样 的 设计 手段 对 飞行 器 的 设计 不 但 效率 极 低 , 而 且 准 确 性 不 够 高 ， 
容易 出 现 飞 行事 故 危及 飞行 员 安 全 ， 甚 至 导致 飞行 器 设计 失败 。 因 此 飞行 器 的 工程 设计 
中 急需 一 种 准确 高 效 的 方法 来 分 析 弹 性 飞行 器 的 跨 音速 气动 弹性 特性 , 并 提供 有 效 克 服 
气动 弹性 问题 的 方法 。 


1.3 本 文 的 主要 工作 


本 文 基于 RANS 方 程 ， 选 取 合 适 的 淇 流 模型 和 算法 ， 定 义 专门 的 边界 条 件 ， 研 制 了 相 
应 的 求解 器 。 具 体 的 工作 内 容 包括 : 
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1. 介绍 了 三 维 直角 坐标 系 下 的 流动 控制 方程 ， 定 义 专门 的 初始 边界 模型 条 件 ， 给 
出 了 基于 有 限 体积 法 下 的 离散 方程 ; 

2. 各 种 注 流 模型 的 对 比 ， 采 用 -a dioi RAT REM, 

3. 详细 介绍 了 三 维 SIMPLE 算 法 ， 给 出 求解 器 的 设计 方法 ， 研 制 了 相应 的 软件 。 
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第 二 章 ”计算 流体 力学 


2.1 CFD 方法 综述 


CFD 方 法 是 对 流 场 的 控制 方程 组 用 数值 方法 将 其 离散 到 一 系列 网 格 节点 上 ， 并 求 其 
离散 数值 解 的 一 种 方法 。 控 制 所 有 流体 流动 的 基本 规律 是 ， 质量 守恒 定律 、 动 量 守恒 定 
律 和 能 量 守 恒定 律 。 由 它们 可 以 分 别 导 出 连续 性 方程 、 动 量 方程 和 能 量 方程 ， 得 到 纳 维 
斯 托 克 斯 偏 微分 方程 组 ,简称 为 N 一 8 方程 组 。N 一 S 方 程 组 是 流体 流动 所 必需 遵守 
的 普遍 规律 。 在 守恒 方程 组 的 基础 上 ， 加 上 反映 流体 流动 特殊 性 质 的 数学 模型 (如 满 流 
模型 、 人 燃烧 模型 、 多 相 流 模型 等 ) 和 边界 条 件 、 初 始 条 件 ， 构 成 封闭 的 方程 组 来 数学 描 
述 特定 流 场 、 流 体 的 流动 规律 。 





尔 


2.2 ”CFD 的 求解 步骤 外 


1 建立 反映 工程 问题 或 物理 问题 本 质 的 数学 模型 。 具 体 地 说 就 是 要 建立 反映 问题 
各 个 量 之 间 关 系 的 微分 方程 及 相应 的 定 解 条 件 ， 这 是 数值 模拟 的 出 发 点 。 没 有 正确 完善 
的 数学 模型 ， 数 值 模拟 就 毫 无 意义 。 流 体 的 基本 控制 方程 通常 包括 质量 守恒 方程 、 动 量 
守恒 方程 、 能 量 守恒 方程 ， 以 及 这 些 方 程 相应 的 定 解 条 件 。 

2 寻求 高 效率 、 高 准确 度 的 计算 方法 ， 即 建立 针对 控制 方程 的 数值 离散 化 方法 ， 
如 有 限 差分 法 、 有 限 元 法 、 有 限 体积 法 等 。 

3 ”编制 程序 和 进行 计算 。 这 部 分 工作 包括 计算 网 格 划分 、 初 始 条 件 和 边界 条 件 的 
输入 、 控 制 参数 的 设 定 等 。 这 是 整个 工作 中 花 时 间 最 多 的 部 分 。 由 于 求解 比较 复杂 ， 比 
如 Navier-Stokes 方程 就 是 一 个 十 分 复杂 的 非 线性 方程 , 数值 求解 方法 在 理论 上 不 是 绝 
对 完善 的 , 所 以 需要 通过 实验 加 以 验证 。 正 是 从 这 个 意义 上 讲 ， 数 值 模 拟 又 叫 数值 试验 。 
应 该 指出 ， 这 部 分 工作 不 是 轻而易举 就 可 以 完成 的 。 

4 显示 计算 结果 。 计 算 结 果 一 般 通 过 图 表 等 方式 显示 ， 这 对 检查 和 判断 分 析 质 量 
和 结果 有 重要 参考 意义 。 
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2. 3 CFD 的 求解 过 程 







建立 控制 方程 
确立 初始 条 件 及 边界 条 件 


划分 计算 网 格 、 生 成 计算 节点 


建立 离散 方程 







离散 初始 条 件 和 边界 条 件 







给 定 求解 控制 参数 







求解 高 散 方 程 









是 

图 2.1 CD 工作 流程 图 

总 体 计算 流程 如 图 2. 1 Bras 

1 建立 控制 方程 

建立 控制 方程 ， 是 求解 任何 问题 前 都 必须 首先 进行 的 。 一 般 来 讲 ， 这 一 步 是 比较 简 
A. 

2 确定 边界 条 件 与 初始 条 件 

初始 条 件 与 边界 条 件 是 控制 方程 有 确定 解 的 前 提 ， 控 制 方程 与 相应 的 初始 条 件 、 边 
界 条 件 的 组 合 构成 对 一 个 物理 过 程 完整 的 数学 描述 。 

初始 条 件 是 所 研究 对 象 在 过 程 开始 时 刻 各 个 求解 变量 的 空间 分 布 情况 。 对 于 瞬 态 问 
题 ， 必 须 给 定 初始 条 件 ， 对 于 稳 态 问题 ， 不 需要 初始 条 件 。 

边界 条 件 是 在 求解 区 域 的 边界 上 所 求解 的 变量 或 其 导数 随地 点 和 时 间 的 变化 规律 。 
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对 于 任何 问题 ， 都 需要 给 定 边界 条 件 。 

3 划分 计算 网 格 

采用 数值 方法 求解 控制 方程 时 ， 都 是 想 办 法 将 控制 方程 在 空间 区 域 上 进行 离散 ， 然 
后 求解 得 到 的 离散 方程 组 。 要 想 在 空间 域 上 离散 控制 方程 ， 必 须 使 用 网 格 。 现 已 发 展 出 
多 种 对 各 种 区 域 进行 离散 以 生成 网 格 的 方法 ， 统 称 为 网 格 生成 技术 。 

”不 同 的 问题 采用 不 同 数值 解法 时 ， 所 需要 的 网 格 形式 是 有 一 定 区 别 的 ， 但 生成 网 格 
的 方法 基本 上 是 一 致 的 。 目 前 ， 网 格 分 为 结构 网 格 和 非 结构 网 格 两 大 类 。 

4 建立 离散 方程 

对 于 在 求解 域内 所 建立 的 偏 微分 ， 理 论 上 是 有 真 解 (或 称 精确 解 或 解析 解 ) 的。 但 
由 于 所 处 理 问题 自身 的 复杂 性 ， 一 般 很 难 获得 方程 的 真 解 。 因 此 ， 就 需要 通过 数值 方法 
把 计算 域内 有 限 数量 位 置 (网 格 节点 或 网 格 中 心 点 ) 上 的 因 变 量 值 当 作 基 本 未 知 量 来 处 
理 ， 从 而 建立 一 组 关于 这 些 未 知 量 的 代数 方程 组 ， 然 后 通过 求解 代数 方程 组 来 得 到 这 些 
节点 值 ， 而 计算 域内 其 他 位 置 上 的 值 则 根据 节点 位 置 上 的 值 来 确定 。 | 

由 于 所 引入 的 应 变量 在 节点 之 间 的 分 布 假设 及 推导 离散 化 方程 的 方法 不 同 , 就 形成 
了 有 限 差分 法 、 有 限 元 法 、 有 限 体积 法 等 不 同类 型 的 离散 化 方法 。 

5 离散 初始 条 件 和 边界 条 件 

前 面 所 给 定 的 初始 条 件 和 边界 条 件 是 连续 性 的 ， 如 在 静止 壁面 上 速度 为 0， 现 在 需 
要 针对 所 生成 的 网 格 ,将 连续 性 的 初始 条 件 和 边界 条 件 转化 为 特定 节点 上 的 值 ， 如 静止 
壁面 上 共有 90 个 节点 ， 则 这 些 节点 上 的 速度 值 应 均 设 为 0。 

在 商用 CFD 软件 中 ， 往往 在 前 处 理 阶段 完成 了 网 格 划分 后 ， 直 接 在 边界 上 指定 初始 
条 件 和 边界 条 件 ， 然 后 由 前 处 理 软件 自动 将 这 些 初始 条 件 和 边界 条 件 按 离 艇 的 方式 分 配 
到 相应 的 节点 上 去 。 

6 给 定 求解 控制 参数 

在 离散 空间 上 建立 了 离散 化 的 代数 方程 组 ， 并 施加 离散 化 的 初始 条 件 和 边界 条 件 
后 ， 还 需要 给 定 流体 的 物理 参数 和 汕 流 模型 的 经 验 系数 。 此 外 ， 还 要 给 定 和 迭代 计算 的 控 
制 精度 、 瞬 态 问题 的 时 间 步 长 和 输出 频率 等 。 

在 CFD 的 理论 中 ,这 些 参数 并 不 值得 去 探讨 和 研究 ， 但 在 实际 计算 时 ， 它 们 对 计算 
的 精度 和 效率 有 着 重要 的 影响 。 
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7 求解 离散 方程 

在 进行 了 上 述 设置 后 ， 生 成 了 具有 定 解 条 件 的 代数 方程 组 。 对 于 这 些 方程 组 ， 数 学 
上 已 有 了 相应 的 解法 ,如 线形 方程 组 可 采用 Gauss 消去 法 或 Gauss-Seidel 迭代 法 求解 ， 
而 对 非 线性 方程 组 ， 可 采用 Newton-Raphson 方法 。 

8 判断 解 的 收 全 性 

对 于 稳 态 问题 的 解 ， 或 是 瞬 态 问题 在 某 个 特定 时 间 步 上 的 解 ， 往 往 需要 通过 多 次 选 
代 才 能 得 到 。 有 时 ， 因 网 格 形式 或 网 格 大 小 、 对 流 项 的 离散 插值 格式 等 原因 ， 可 能 导致 
解 的 发 散 。 对 于 咀 态 问题 ， 若 采用 显 式 格 式 进行 时 间 域 上 的 积分 ， 当 时 间 步 长 过 大 时 ， 
也 可 能 造成 解 的 振 功 或 发 散 。 因 此 ， 在 迭代 过 程 中 ， 要 对 解 的 收敛 性 随时 进行 监视 ， 并 
在 系统 达到 指定 精度 后 ， 结 束 迭 代 过 程 。 

9 显示 和 输出 计算 结果 

通过 上 述 求解 过 程 得 出 了 各 计算 点 上 的 解 后 , 需要 通过 适当 的 手段 将 整个 计算 域 上 
的 结果 表示 出 来 。 这 时 ， 我 们 可 采用 线 值 图 、 矢 量 图 、 等 值 线 图 、 流 线 图 、 云 图 等 方法 
对 计算 结果 进行 表示 
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第 三 章 流体 运动 的 控制 方程 和 泣 流 模型 


流体 运动 现象 大 量 的 出 现在 自然 界 及 各 个 工程 领域 中 ， 其 具体 的 表现 形式 多 种 多 
样 。 任何 流体 运动 的 动力 学 特性 都 受到 最 基本 的 三 个 定律 ， 即 质量 守恒 定律 、 动 量 守 恒 
定律 和 能 晤 守恒 定律 的 支配 。 质 量 守恒 定律 :单位 时 间 内 微 元 体 中 流体 质量 的 增加 ， 等 
于 同一 时 间 间 隔 内 流入 该 微 元 体 的 净 质 量 。 动 量 守 恒定 律 : 微 元 体 中 流体 动量 的 增加 率 
等 于 外 界 作用 在 该 微 元 体 上 的 各 种 力 之 和 。 能 量 守恒 定律 : 微 元 体 中 能 量 的 增加 率 等 于 
进入 微 元 体 的 净 热 流量 加 上 体积 力 和 表面 力 对 微 元 体 所 做 的 功 。 这 些 基 本 定律 可 由 数学 
方程 组 来 决定 。 根 据 上 述 三 个 定律 ， 可 以 得 到 飞行 器 周围 流体 运动 的 控制 微分 方程 的 矢 
量 形式 为 & 7, 8, 9, ". 


0 
1 连续 方程 : :+ Yi(OF) = 0 (3.1) 
2. zjB7; Pe: Oe ey ye ye np oe (3.2) 
of & y & pu 
. (0l) oy 
3. 能 量 方程 : —— +V- (piu)  -pV-u e V (kVT)4 95, (3.3) 


at 

其 中 jp 为 流体 密度 ,u 为 流体 速度 ，k 为 流体 热传导 率 ，% RARAN, i 是 单位 质量 气 
TERT. FE RTS ERR AEF ons 都 等 于 0 

流体 力学 中 有 五 个 基本 参量 ;三 个 速度 分 量 和 两 个 热力 学 参量 。 例 如 压强 、 温 度 、 
密度 、 答 和 焙 中 的 任何 两 个 热力 学 参量 ， 都 足以 确定 流体 的 热力 学 状态 ， 因 此 也 可 以 确 
定 出 所 有 其 他 的 热力 学 参量 。 一 旦 规定 了 速度 矢量 V 和 两 个 热力 学 参量 随时 间 和 空间 变 
化 的 函数 ， 我 们 就 可 以 完全 确定 流体 的 流 场 。 因 此 ， 我 们 需要 五 个 独立 的 方程 。 通 常 是 
三 个 运动 方程 的 分 量 方程 , 一 个 连续 方程 和 一 个 能 量 方程 。 为 了 便于 用 三 个 变量 (温度 、 
密度 和 压强 ) 来 写 出 能 量 方程 往往 还 引入 一 个 物 态 方程 。 在 这 种 情形 下 ， 一 共有 六 个 
变量 和 六 个 方程 。 从 而 得 出 三 维 定常 可 压缩 流体 守恒 型 数学 模型 ”: 
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ap la) al) (ow) wop 
oy = ny — . = 

连续 roar? = UY (pV)=0 (3. 4) 
X 动 量 ; ae AED, div(pu) = - divin gradu) - P +, (3.5) 
Y 动量 : Maur aot] Bot (3.6) 
Z 动量 : Nm, div(pwV) = div(u grad v)- Tus, (3. 7) 

alei) | s.s as V. 
能 量 方程 tv (piu) = —pV-u+V-(kVT)+9+5, (3.8) 
RAGE: P-P(oT)  i=i(p,T) (3.9) 


(完全 气体 状态 方程 为 ，p = pRT,i=c,T ) 
方程 组 数目 和 未 知 函 数 数目 相等 ， 在 数学 上 这 个 函数 是 封闭 的 ， 即 如 果 提 供 适 
当 的 辅助 条 件 这 个 方程 组 是 可 解 的 。 
为 方便 看 到 各 个 守恒 方程 的 共同 点 引入 流体 物理 量 的 一 般 表 达 式 g ,所 有 守恒 形 
式 方程 可 写 为 如 下 有 用 形式 ， 即 : 
“08 + div( pug) = div(T grad 4) S (3.10) 
ERKKA: HAM (transient term), Wi (convective term), 7 BAH 
(diffusive term), W (source term)， 用 文字 表述 为 : 流体 元 的 4 的 增加 率 加 
上 $$ 的 净 流 出 率 等 于 gp 的 扩散 增加 率 加 上 y 的 源 增 加 率 ， 写 出 求解 方程 (3. 10) 的 源 
程序 ， 对 于 不 同 的 pg， 只 要 重复 调用 该 程序 ， 并 给 定 TT 和 S$ 的 适当 表达 式 以 及 适当 的 
初始 条 件 和 边界 条 件 ， 即 可 求解 。 


3. 2 ” 满 流 对 时 均 NS 方程 的 影响 


在 自然 界 和 工程 实际 中 , 最 经 常 发 生 的 流动 状态 是 江 流 , 这 时 的 流动 不 再 是 平滑 的 ， 
而 始终 存在 着 很 不 规则 的 脉动 。 满 流 流动 随时 间 和 空间 都 呈现 出 不 规则 的 脉动 ， 它 们 是 
由 大 大 小 小 的 游 涡 引起 的 。 淇 流 是 紊乱 的 ， 这 是 淇 流 的 本 质 。 淇 流 不 是 完全 随机 的 ， 因 
为 潮流 也 服从 流体 运动 基本 方程 组 ， 如 果 假 设 某 一 速度 分 量 是 完全 随机 的 ， 则 其 余 两 个 
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分 量 一 定 会 由 三 大 守恒 定律 限制 其 脉动 的 范围 。 研 究 发 现 ， 在 湛 流 混合 层 和 边界 层 中 都 
存在 着 拟 序 结构 ， 它 们 都 以 大 尺度 尝 涡 流动 为 特征 。 流体 运动 作为 由 Navier-Stokes Jj 
程 组 描写 的 确定 性 过 程 ， 在 高 雷诺 数 时 会 出 现 类 似 于 随机 的 ， 非 确定 性 的 物理 特性 。 

工程 中 过 到 的 流动 大 多 数 都 是 浇 流 ， 通 常 认为 ， 无 论 江 流 运动 多 么 复杂 ， 非 稳 态 的 
连续 方程 和 Navier-Stokes 方程 仍然 适用 于 灌流 的 有 瞬时 运动 。 

灌流 运动 看 作 由 时 间 平 均 流 动 和 有 瞬时 脉动 流动 全 加 而 成 。 

任意 变量 gp 的 时 间 平 均值 定义 为 ”: 


一 1 +åt 
Pf g(t yt (3.11) 


这 里 ， 上 标 “ 一 ”代表 对 时 间 的 平均 值 。 如 果 用 上 标 “ ”代表 脉动 值 ， 物 理 量 的 
瞬时 值 pg 、 时 均值 g 及 脉动 值 多 之 间 有 如 下 关系 ; 





$-6 4 (3. 12) 
现在 用 平均 值 与 脉动 值 之 和 代替 流动 变量 ， 即 
VeV«V; ueiu; VeV+tv; waW+w; pepep (3. 13) 
在 考虑 脉动 的 情况 下 ， 把 (3-13) 代入 (3-4) (3-5) (3-6) (3-7) PB BT 
均 流 动 的 控制 方程 如 下 : 
divū = 0 (3. 14) 
Qu^ Quv Ww 
a P)» earam | -区 x J (3. 15) 
La x EU .9uv av? avw 
Z edi GP) 5E svangan uc a a (3. 16) 
qb) D oras f | (3.17) 
à p à à ay a 
对 于 其 他 变量 % 的 输 运 方程 (如 能 量 方 程 ) 作 类 似 处 理 ， 可 得 
PP + div(GV) =div(T grad )+ |. ne. us s (3.18) 


但 实际 流动 中 ， 密 度 可 能 是 变化 的 。 在 此 ， 忽略 密 度 脉动 的 影响 ,但 考虑 平均 密度 
的 变化 ， 写 出 可 压 游 流 平均 流动 的 控制 方程 如 下 注意， 为 方便 起 见 ， 除 脉动 值 的 时 均 
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值 外 ， 下 式 中 去 掉 了 表示 时 均值 的 上 划 线 符合 “一 ”， 如 # 用 5 来 表示 ) 
连续 方程 ; $E p divov) = =0 (3. 19) 
动量 方程 (Navier-Stokes 方程 ): 


A * div(puV) = div(u grad u)~ Z+ Bue Rue «s ^ (3. 20) 





ax oy az 
2 rt 
SO dispu)» divin grad v) -E+ | xa reni xe. : (3.21) 
n 

UPW), aivgwV) td m- P. ee x 22. S, (3.22) 
其 他 变量 的 输 运 方程 ， 
(pp). .. ;  &(ou'g) opr$) 3pw) 2], 

x: + div(pgV) = div(T grad $) + | s y A (3. 23) 


A (3.190 为 时 均 形 式 的 连续 方程 ， 方 程 (3.20) (3.21) (3.22) 为 时 均 形式 的 
Navier-Stokes 方程 。 由 于 采用 的 是 Reynolds 平均 法 , 因此 , 方程 (3.20)(3.21)(3.22) 
被 称 为 Reynolds 时 均 Navier-Stokes 方程 (Reynolds-Averaged Navier-Stokes， 简 称 
RANS)， 常 直接 称 为 Reynolds 方程 。 可 以 看 到 ， 时 均 流动 的 方程 中 多 出 与 -pr 有 关 的 
项 ， 我 们 定义 该 项 为 Reynolds 应 力 ， 即 : 

Ty = -puu (3.24) 

RB, cQ EBM 6 个 不 同 的 Reynolds MAT, BB 3 个 正 应 力 和 3 个 切 应 力 。 

Fast (3. 19) (3.20) (3.21) (3.22) (3.23) 构成 的 方程 组 共有 五 个 方程 ， 现 在 新 
增 了 式 〈3. 24) , 即 6 个 Reynolds 应 力 ， 再 加 上 原来 的 5 个 时 均 未 知 量 Cus vw, p.62, 
总 共有 11 个 未 知 量 ， 因 此 ， 方 程 组 不 封闭 ， 必 须 引入 新 的 汕 流 模型 方程) 才能 使 方 
程 组 封闭 。 要 使 方程 组 封闭 ， 就 得 把 漠 流 的 脉动 值 与 时 均值 等 联系 起 来 ， 这 就 必须 对 
Reynolds 应 力 做 出 某 种 假定 。 
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3. 3 ”灌流 模型 的 比较 


根据 对 Reynolds 应 力 做 出 的 假定 或 处 理 方式 不 同 , 目前 常用 的 洋流 模型 有 两 大 类 : 
Reynolds 应 力 模型 和 涡 粘 模型 ， 如 图 3. 1 所 示 LL EN 


混合 长 度 模 型 
I. 
一 方程 模型 标准 ~ 模型 


- J 
s ded Realizable* -£ 
* m 
典 模型 ia 


雷诺 应 力 方程 模型 
雷诺 应 力 模型 


代数 应 力 方程 模型 


图 3. 1: 基于 雷诺 方程 的 经 典 模型 结构 示意 图 


(1) Reynolds 应 力 模 型 

在 Reynolds 应 力 模 型 计算 量 很 大 , 这 就 极 大 地 限制 了 Reynolds 应 力 模 型 在 工程 问 
题 中 的 应 用 。 因 此 这 里 不 做 介绍 

(2) 涡 粘 模型 5 

在 涡 粘 模型 方法 中 ， 不 直接 处 理 Reynolds SJA, MESA mA turbulent 
viscosity)， 或 称 涡 粘 系数 (eddy viscosity)， 然 后 把 庙 流 应 力 表示 成 测 动 粘度 的 函 
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数 ， 整 个 计算 的 关键 在 于 确定 这 种 灌 动 粘度 。 
ae HE RE AGE RAF Boussinesq 提出 的 涡 粘 假定 , 该 假定 建立 了 Reynolds WH 
方程 ， 即 : 


p Qu; au, 2 Ou, 
puu’ = —4—]-—] ok + u — Iô. (3. 25) 
ds Zr ox, | | 2 ax, | : 


这 里 ，j Aim, u 为 时 均 速 度 ，6, Rb “Kronecker derta” W5 (4i jf, 
6,=1; Siw jh, 6,=0), k AmaE (turbulent kinetic energy): 


上 -人 (3.26) 
2 “2 


HOWE p, 是 空间 坐标 的 函数 ， 取 决 于 流动 状态 ， 而 不 是 物性 参数 。 由 上 可 见 ， 引 
入 Boussinesq 假定 以 后 ， 计 算 测 流 流动 的 关键 就 在 于 如 何 确定 上 。 这 里 所 请 的 涡 粘 模 
型 ， 就 是 把 与 清流 时 均 参 数 联系 起 来 的 关系 式 。 

依据 确定 u 的 微分 方程 数目 的 多 少 ， 涡 粘 模型 包括 ， 零 方程 模型 、 一 方程 模型 、 两 
方程 模型 。 


3.3.1 零 方 程 模型 


所 渭 零 方程 模型 是 指 不 使 用 微分 方程 ， 而 是 用 代数 关系 式 ， 把 潮 动 粘度 与 时 均值 联 
系 起 来 的 模型 。 它 只 用 满 流 的 时 均 连 续 方程 ， Reynolds 方程 组 成 方程 组 ， 把 方程 组 中 的 
Reynolds 应 力 用 平均 速度 场 的 局 部 速度 梯度 来 表示 。 
(1) 混合 长 度 模 型 (mixing length model) 
零 方 程 模型 方案 有 多 种 ， 最 著名 的 是 Prandtl 提出 的 混合 长 度 模型 (mixing length 
model). Prandtl 假定 渤 动 粘度 u, 正比 于 时 均 速 度 志 的 梯度 和 混 介 长 度 RRI. 
混合 长 度 理论 的 优点 是 直观 简单 ， 对 于 如 射流 、 混 合 层 、 扰 动 和 边界 层 等 带 有 沙 的 
前 切 层 的 流动 比较 有 效 , 但 只 有 在 简单 流动 中 才 比 较 容易 给 定 混合 长 度 1, ， 对 复杂 流动 
则 很 难 确定 1。， 而 且 不 能 用 于 模拟 带 有 分 离 及 回流 的 流动 ， 因 此 ， 零 方程 模型 在 实际 工 
程 中 很 少 使 用 。 
(2) Baldwin-lomax (BL) Bi78 P-99 
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在 复杂 的 流动 数值 模拟 (N-S 方程 ) 中 所 使 用 的 皆 是 简化 的 满 流 模型 ， 尤 以 零 方程 一 
一 BL 模型 使 用 更 为 广泛 
BL 模型 较 适 合 于 平衡 状态 的 油 流 流动 , 包括 附着 流动 , 或 具有 较 小 分 离 气泡 的 分 离 
流动 。Degani 等 在 研究 三 维 炮 弹 型 物体 的 超声 速 潮流 流动 中 发 现 ， 采 用 BL 模型 ， 在 小 
迎 角 (附着 流 ) 时 ,计算 结果 与 实验 结果 很 一 致 。 但 在 大 迎 角 下 ( 即 存在 横向 分 离 时 ) ， 两 
者 不 一 致 。 


3.3.2 一 方程 模型 


在 零 方 程 模型 中 ， 淇 动 粘度 j MIA KIEL, 都 把 Reynolds 应 力 和 当地 平均 速度 梯 
度 相 联 系 ， 是 一 种 局 部 平衡 的 概念 ， 忽 略 了 对 流 和 扩散 的 影响 。 为 了 弥补 混合 长 度假 定 
的 局 限 性 ， 人 们 建议 在 满 流 的 时 均 连 续 方程 和 Reynolds 方程 的 基础 上 ， 再 建立 一 个 注 
动能 上 的 输 运 方程 ， 而 表示 成 x 的 函数 ， 从 而 可 使 方程 组 封闭 。 

这 里 ， 沸 动能 k MMAR SAM: 


3 
alok) aloku) a f(a, yp, ( ,oo 
à ^ ae, — de Ges "A ax, ar, 3; (9, oD 


i i i 








MEBA, HE ESKUCAEESIR. SRT. P RON. PAT. BBO. m 

Kolmogorov-Prandtl 表达 式 有 : 
I = pC, vkl (3. 28) 

EP, oQCQC, 为 经 验 常数 ， 多 数 文献 建议 os -LC, ~0.09, MC, 的 取 值 在 不 
同 的 文献 中 结果 不 同 ， 从 0. 08 到 0. 38 不 等 。 但 这 个 问题 在 后 面 要 介绍 的 两 方程 模型 中 
不 存在 。! 为 潮流 脉动 的 长 度 比 尽 ， 依 据 经 验 公 式 或 实验 而 定 o 

X (3.27) 和 (3.28) 构成 一 方程 模型 。 一 方程 模型 考虑 到 满 动 的 对 流 输 运 和 扩散 
输 运 ， 因 而 比 零 方 程 模型 更 为 合理 。 但 是 一 方程 模型 的 缺点 是 如 何 确定 长 度 比 尺 1 仍 为 
不 易 解决 的 问题 ， 因 此 很 难得 到 推广 应 用 。 


3. 3.3 ”两 方程 模型 
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在 两 方程 模型 中 ， 比 较 典 型 的 是 标准 k -e 模型 ， 是 在 一 方程 模型 的 基础 上 ， 新 引入 
一 个 关于 汝 流 耗 散 率 e 的 方程 后 形成 的 。 大 -sg 模型 是 针对 滑 流 发 展 非常 充分 的 湛 流 流动 
来 建立 的 ， 也 就 是 说 ， 它 是 一 种 针对 高 Re Sm VLL ARH, Re 数 较 低 的 流动 使 用 
k-e 模型 进行 计算 ， 就 会 出 现 问题 。 这 时 ， 必 须 采 用 特殊 的 处 理 方 式 ， 来 解决 近 壁 区 
的 流动 计算 及 低 Re 数 时 的 流动 计算 问题 。 常 用 的 解决 方法 有 两 种 ， 一 种 是 采用 壁面 函 
数 法 ， 另 一 种 是 采用 低 Re 数 的 ke 模型 。 标 准 t-e 模 型 比 零 方程 模型 和 一 方程 模型 有 了 
很 大 改进 , 在 科学 研究 及 工程 实际 中 得 到 了 最 为 广泛 的 检验 和 成 功 应 用 , 但 用 于 强 旋 流 、 
弯曲 壁面 流动 或 弯曲 流 线 流 动 时 ， 会 产生 一 定 的 失真 。 该 模型 是 目前 使 用 最 广泛 的 测 流 
模型 。 

ERIE, AEk -e 模型 在 计算 高 Re 数 满 流 方 面具 有 广泛 的 应 用 , 因此 本 文采 取 
标准 上 -< 模型 作为 潮流 模型 
3. 4 k-e 模型 

ERT AR k 的 方程 的 基础 上 ， 再 引入 一 个 关于 潢 动 耗 散 率 e 的 方程 ， 便 形成 了 
k-e 两 方程 模型 ， 称 为 标准 x- 模 型， 在 模型 中 ， 表 示 漳 动 耗 散 率 的 被 定义 为 : 
EE (3. 29) 
MOE uw, 可 表示 成 -zs 的 函数 ， 即 

H, - pc, E. (3. 30) 
€ 

对 可 压缩 流体 采用 AEk -e 模型 ， 方 程 和 方程 分 别 表示 为 :; 


a à ðu , 
alok) , aloku) 3 u+ x +H; kl Baa d wo +G, -pe-Y, +S, (3.31) 
ar à, ax, c, ) a, ax, a, Jax, 











a(p yy (peu) a poe ae +H, IE 2 ec, (6 «C,6,) 
at à — dx e, Jàx, Ox, ji, k 


(3. 32) 


2 

€ 
-C,p-—-*$ 
2, P k E 


其 中 ，G, RAP PRR SoH k 的 产生 项 ， G, ECT VEU S EG RO 
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动能 k 的 产生 项 ，J fo HD Hi Vien BK SKI VERA, CC, ,C3, 为 经 验 常 数 , ok,ar。 
分 别 是 与 潮 动 能 K 和 耗 散 率 对 应 的 Prandt1 数 ，8 ,9 是 用 户 定义 的 源 项 , 当 不 考虑 自 
定义 源 项 时 ，S =5, =0。 
各 项 的 计算 公式 如 下 中 ， 
首先 ，G, 是 由 于 平均 速度 梯度 引起 的 满 动能 k 的 产生 项 ， 由 下 式 计算 : 


G, LEES (3.33) 


式 (3.33) 的 展开 式 为 


2 2 2 2 2 2 
sup) (2) (2) HE) (E e | 03.34) 
ax oy 0z By dx oz Ox oz y 

G, 是 由 于 浮力 引起 的 满 动 能 k 产生 项 ， 对 于 可 压 流体 ， 有 ; 


4, oT 
G, = Bg, —— (3. 35) 
>= PB Pr, ox, 


其 中 , Pr deta) Prandtl 数 ， 在 该 模型 中 可 取 Pr =0.85, g, 是 重力 加 速度 在 第 i y 
向 的 分 量 ， 是 热 脱 胀 系数 ， 可 由 可 压 流体 的 状态 方程 求 出 ， 其 定义 为 : 


= (3. 36) 
p aT 
Y, 代表 可 压 满 流 中 脉动 扩张 的 贡献 ， 对 于 可 压 流 体 ， 有 : 
Y, = 2peM,’ (3. 37) 


其 中 ，M， 是 消 动 Mach 数 ，M, = £, 4 是 声速 ，a = VRT 。 


式 (3-31)，《〈3-32) 中 的 各 系数 为 经 验 常数 ， 取 值 如 表 3. 1: 
表 3.1 经 验 常数 


当主 流 方向 与 重力 方向 平行 时 ， cx =1， 当 主流 方向 与 重力 方向 垂直 时 ，c,, =0 
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3. 5 小 结 


本 章 主 要 给 出 了 三 维 直角 坐标 下 的 流 场 控制 方程 ， 比 较 了 不 同 的 测 流 模型 的 优 缺 
Ao BAART E-E BULUM, IB T ke 数学 模型 推导 和 经 验 系数 
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第 四 章 ”边界 条 件 


在 求解 RANS 方程 组 时 ， 要 想 方 程 组 的 解 具有 唯一 性 ， 必 须要 有 一 定 的 初始 条 件 和 
边界 条 件 ， 因 此 ， 确 定 的 边界 条 件 是 RANS 方程 组 数值 求解 的 重要 问题 之 一 。 边 界 条 件 
的 处 理 方法 在 很 大 程度 上 影响 着 数值 计算 过 程 中 的 稳定 性 ， 如 果 处 理 不 当 ， 往往 导 致 数 
值 计 算 不 能 收敛 , 甚至 发 散 。 另 外 , 边界 条 件 的 处 理 往往 对 流 场 精度 有 非常 重要 的 影响 。 
为 有 效 提 高 流 场 计算 的 精度 ， 有 必要 提高 边界 格式 的 处 理 精度 ， 使 格式 在 边界 上 和 场 内 
精度 相 匹 配 。 边 界 条件 主 要 包括 流动 进口 边界 ， 流 动 出 口 边 界 ， 壁 面 边 界 ， 对 称 边 界 。 


4.1 iR& mg 


对 于 粘性 绕 流 , 在 物 面 处 , 数值 计算 中 使 用 无 滑 移 、 绝热 壁 及 法 向 零 压力 梯度 条 件 * 


19,20, 21, 2247), 

u=vaweak =Q (4.1) 

其 中 心 ，v，m 为 物 面 运动 速度 ， 在 固 壁 上 不 能 显 式 给 出 耗 散 率 e 的 条 件 ， 因 为 在 
壁面 上 它 并 不 为 等 。 

考虑 到 理想 气体 状态 方程 :P = pRT , 压力 和 温度 均 满足 法 向 零 梯度 的 条 件 , 因此 ， 


密度 在 物 面 也 满足 法 向 零 梯度 条 件 : “2-0 ， 壁面 附近 的 网 格 比较 密 ， 认 为 附 面 层 
RERE: m 一 0 
on 


4.2 ”其 它 边 界 条 件 


在 有 限 域内 求解 复杂 外 形 绕 流 问题 时 ， 除 了 以 上 壁面 边界 条 件 外 ， 一般 还 涉及 以 下 
几 类 处 理 方式 相近 的 边界 条 件 : 

飞行 器 在 不 考虑 发 动 喷 流 的 条 件 下 ， 常 用 的 边界 条 件 有 ， 

(D 进口 边界 : 进口 边界 上 固定 为 自由 来 流 参数 ，u v w，p 和 TT 是 预先 给 定 的 ， 
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流动 出 口 边界 -一般 选 在 离 几 何 扰动 足够 远 的 地 方 来 施加 ， 也 就 是 说 ， 沿 流动 方向 个 所 有 
变量 4〈 除 压力 外 ) 是 无 变化 的 ， 即 ~0， 消 动能 k RRDKCICROR e 也 是 预先 给 定 的 ， 
满 动能 k 和 汕 流 耗 艇 率 e :一 般 根据 实验 数据 得 出 ， 或 者 是 通过 估算 公式 给 出 。 对 于 平 
fbit: 
k, -S (Us) ieee, 1-007L (4.2) 
sth, 元 是 进口 的 平均 速度 , OK L EORUM, T NNAIR 

(2) 出 口 边界 ， 即 远 场 边界 ， 远 场 理论 上 无 限 远 ， 未 受到 任何 扰动 ， 实 际 上 ， 远 
场 取 有 限 大 的 区 域 ， 以 满足 工程 计算 需要 。 

(3) 对 称 边界 : 对 于 只 要 求 计算 飞行 器 纵向 空气 动力 特性 的 情形 ， 通 常 沿 飞行 器 
纵向 对 称 面 取 半 模 来 进行 计算 分 析 , 这 时 飞行 器 的 纵向 面 就 按 对 称 面 边界 处 理 , 一 般 地 ， 
若非 特殊 说 明 ， 文 中 计算 域 关于 x-y 平面 对 称 。 对 于 全 飞行 器 ， 不 予 考虑 。 对 称 面 边界 
上 的 法 向 速度 为 零 , w = 0， 沿 着 对 称 面 法 向 的 任意 标量 物理 量 的 梯度 为 零 ， 即 对 称 


面 上 通 量 为 零 : 凶 -0 
on 


21 


中 北大 学 学 位 论文 
第 五 章 RANS 方程 离散 化 


5. 1 离散 方法 综述 


由 于 RANS 方程 是 非 线性 的 二 阶 偏 微分 方程 组 ， 所 以 ， 直 接 对 该 方程 求解 析 解 是 
非常 困难 的 。 只 有 极 少数 边界 条 件 简 单 的 流动 ， 方 程 通过 简化 之 后 ， 才 有 可 能 得 到 解析 
解 。 迄 今 为 止 ， 在 实际 工程 问题 中 ， 一 般 是 对 RANS 方程 进行 区 域 离散 化 ， 变 成 一 组 线 
性 代数 方程 ， 然 后 在 计算 机 上 计算 得 到 问题 的 数值 解 。 

所 谓 区 域 离散 化 实质 上 是 用 一 组 有 限 个 离散 点 来 代替 原始 的 连续 空间 。RANS 方 
程 的 数值 离散 方法 主要 有 已 2426: 有 限 体积 法 (Finite Volume Method, FVYM) 有 限 差分 法 
(Finite DifferenceMethod, FDM )、 有 限 元 法 (Finite Element Method, FEM )。 

(1) 有 限 差分 法 

有 限 差分 法 (Finite Difference Method， 简 称 FDM) 是 数值 解法 中 最 经 典 的 方法 。 
它 是 将 求解 域 划分 为 差分 网 格 ， 用 有 限 个 网 格 节点 代替 连续 的 求解 域 ， 然 后 将 偏 微分 方 
程 〈 控 制 方程 ) 的 导数 用 差 商 代替 ， 推 导出 含有 离散 点 上 有 限 个 未 知 数 的 差分 方程 组 。 
求 差 分 方程 组 〈 代 数 方程 组 ) 的 解 ， 就 是 微分 方程 定 解 问题 的 数值 近似 解 ， 这 是 一 种 直 
接 将 微分 问题 变 为 代数 问题 的 近似 数值 解法 。 

FDM 建立 在 经 典 的 数学 逼近 理论 的 基础 上 ， 简 单 且 易 为 人 们 所 接受 ， 处 理 效率 较 
高 。 另 一 方面 ， 控 制 方程 表达 了 质量 守恒 的 物理 定律 ， 而 差分 方程 有 时 不 能 严格 保持 守 
恒 性 质 ， 数 值 解 会 出 现 质量 、 动 量 不 平衡 的 守恒 误差 。 经 典 的 FOM 常常 也 不 能 用 来 正 
确 计 算 间 断 解 。 由 于 使 用 泰勒 级 数 展开 ，FDM 一 般 只 能 用 于 处 理 规则 的 几何 区 域 和 正 
HAMA. FOM 通常 在 数值 解 精度 和 处 理 不 规则 物理 域 方 面 ， 存 在 着 根本 性 的 困难 。 
OD 有限 元 法 

有 限 元 法 (Finite Elenent Method， 简 称 FEM) 与 有 限 差分 法 都 是 广泛 应 用 的 流体 
动力 学 数值 计算 方法 。 有 限 元 法 是 将 一 个 连续 的 求解 域 任意 分 成 适当 形状 的 许多 微小 单 
元 ， 并 于 各 小 单元 分 片 构造 插值 函数 ， 然 后 根据 极 值 原理 ( 变 分 或 加 权 余 量 法 )， 将 问 
题 的 控制 方程 转化 为 所 有 单元 上 的 有 限 元 方程 ， 把 总 体 的 极 值 作 为 各 单元 极 值 之 和 ， 即 


22 


中 北大 学 学 位 论文 


将 局 部 单元 总 体 合 成 ， 形 成 嵌入 了 指定 边界 条 件 的 代数 方程 组 ， 求 解 该 方程 组 就 得 到 各 
节点 上 待 求 的 函数 值 。 有 限 元 法 的 基础 是 极 值 原理 和 划分 插值 ， 它 吸收 了 有 限 差分 法 中 
离散 处 理 的 内 核 ， 又 采用 了 变 分 计算 中 选择 逼近 函数 并 对 区 域 进行 积分 的 合理 方法 ， 是 
这 两 类 方法 互相 结合 、 取 长 补 短发 展 的 结果 。 它 具有 很 广泛 的 适应 性 ， 特 别 适 用 于 几何 
及 物理 条 件 比 较 复杂 的 问题 ， 而 且 便 于 程序 的 标准 化 。 对 椭圆 型 方程 有 更 好 的 适用 性 。 
有 限 元 法 因 求解 速度 较 有 限 差分 法 和 有 限 体积 法 慢 ,在 运行 FEM 进行 CFD 数值 计算 时 ， 
每 一 步 计 算 都 需要 求解 一 个 大 型 线性 带 状 三 对 角 方 程 组 ， 需 要 耗费 大 量 机 时 ， 即 使 该 方 
程 组 是 显 式 的 ， 求 解 效 率 也 不 高 。 

(3) 有 限 体积 法 

， 有限 体积 法 (Finite Volume Method， 简 称 FVM)， 是 近年 来 发 展 非常 迅速 的 一 种 离 
散 化 方法 ， 其 特点 是 近似 效率 高 。 

有 限 体积 法 (Finite Volume Method〉 又 称 为 控制 体积 法 (Control Volume Method, 
简称 CVM)。 其 基本 思路 是 : 将 计算 区 域 划分 为 网 格 ， 并 使 每 个 网 格 点 周围 有 一 个 互 不 
重复 的 控制 体积 : 将 待 解 微分 方程 〈 控 制 方程 》 对 每 一 个 控制 体积 积分 ， 从 而 得 出 一 组 
离散 方程 。 其 中 的 未 知 数 是 网 格 点 上 的 因 变 量 几 。 为 了 求 出 控制 体积 的 积分 ， 必 须 假定 
o 值 在 网 格 点 之 间 的 变化 规律 。 从 积分 区 域 的 选取 方法 看 来 ， 有 限 体积 法 属于 加 权 余 量 
法 中 的 子 域 法 ， 从 未 知 解 的 近似 方法 看 来 ， 有 限 体 积 法 属于 采用 局 部 近似 的 离散 方法 。 
简 言 之 ， 子 域 法 加 离散 ， 就 是 有 限 体积 法 的 基本 方法 。 

FVM 类 似 FEM， 是 将 计算 域 划分 成 若干 规则 或 者 不 规则 形状 的 单元 或 者 控制 体 。 
在 计算 出 通过 每 个 控制 体 边 界 沿 法 向 输入 、 输 出 的 质量 通 量 和 动量 通 量 后 ， 对 每 个 控制 
体 分 别 进 行 质量 和 动量 平衡 计算 ， 得 到 线性 方程 组 。 因 此 ，FVM 正 是 对 于 推导 原始 微 
分 方程 所 用 的 控制 体 思想 的 回归 , 与 FDM 和 FEM 的 数值 逼近 相 比 ,其 物理 意义 更 加 直 
接 明晰 。 因 为 两 个 相 邻 控制 体 之 间 的 控制 面 上 输送 的 通 量 ， 对 于 相 邻 控制 体 来 说 大 小 相 
等 、 方 向 相反 ， 因 此 对 于 整个 计算 域 而 言 ， 沿 所 有 内 部 边界 的 通 量 互相 抵消 。 对 于 由 一 
个 或 多 个 控制 体 所 组 成 的 任意 区 域 , 以 及 整个 求解 域 , 都 能 非常 严格 地 满足 物理 守恒 律 ， 
不 存在 守恒 误差 ， 而 且 能 正确 计算 间断 。 

由 于 FVM 继承 了 FDM 的 丰富 格式 ， 具 有 良好 的 守恒 人 性， 能 像 FEM 那样 采用 各 种 
形状 的 网 格 以 适应 复杂 的 边界 几何 形状 ， 却 比 FEM 简便 得 多 ， 目 前 在 CFD 领域 得 到 了 
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广泛 应 用 ， 大 多 数 商 用 CFD 软件 都 采用 这 种 方法 。 


5. 2 ”离散 策略 


鉴于 FVM 的 优势 和 应 用 前 景 ， 本 文采 用 FVM 进行 方程 的 离散 。 


5.3 GT BIRIEFUEISIESDS FER Rn mm 


在 流 场 计算 之 前 ， 首 先 要 将 计算 区 域 进行 离散 化 ， 就 是 对 空间 上 的 计算 区 域 进行 划 
分 ， 把 计算 域 划 分 成 许多 个 子 区 域 ， 并 确定 每 个 区 域 的 节点 ， 生 成 网 格 。 然 后 将 控制 方 
程 在 网 格 上 离散 ， 把 偏向 分 格式 的 控制 方程 转化 为 各 个 节点 上 的 代数 方程 。 

在 求解 域 中 建立 的 偏 微分 方程 ， 理 论 上 是 存在 解析 解 的 。 但 是 由 于 所 处 理 问题 本 身 
的 复杂 性 ， 尤 其 对 于 非 线性 方程 组 ， 很 难得 到 方程 的 解析 解 。 因 此 ， 就 需要 通过 数值 计 
算 的 方法 把 计算 域内 的 网 格 节点 上 的 因 变 量 当 作 基本 未 知 量 来 处 理 ， 从 而 建立 关于 这 些 
未 知 量 的 代数 方程 组 ， 然 后 通过 求解 这 个 代数 方程 组 以 得 
到 节点 上 的 数值 解 。 计 算 域 其 他 位 置 上 的 数值 根据 节点 上 的 值 来 确定 。 

有 限 体 积 法 (Finite Volume Method) 又 称 为 控制 体积 法 (Control Volume 
Method) 。 在 有 限 体积 法 中 将 所 计算 的 区 域 划分 为 一 系列 控制 容积 ， 每 个 控制 容积 都 有 
一 个 节点 作为 代表 ， 通 过 把 控制 方程 对 控制 容积 做 积分 来 导出 离散 力 一 程 。 在 导出 的 离 
散 方 程 中 ,需要 对 界面 上 的 被 求解 函数 本 身 及 其 一 阶 导 数 的 构成 做 出 假设 ， 这 种 构成 的 
方式 就 是 有 限 体积 法 中 的 离散 格式 。 从 积分 区 域 的 选取 方法 来 看 ， 有 限 体积 法 属于 加 权 
余 量 法 中 的 子 域 法 :从 未 知 解 的 近似 方法 看 来 , 有 限 体 积 法 是 属于 局 部 近似 的 离散 方法 。 
所 以 说 子 域 法 加 上 离散 就 是 有 限 体积 法 的 基本 方法 。 因 为 有 限 体 积 法 得 出 的 离散 方程 要 
求 因 变量 的 积分 守恒 对 任意 一 个 控制 体积 都 得 到 满足 ， 对 整个 计算 域 自然 也 得 到 满足 ， 
所 以 采用 有 限 体积 法 的 离散 方程 可 以 保证 具有 守恒 性 ， 这 是 有 限 体积 法 的 优点 。 例 如 有 
限 差 分 法 ， 仅 在 网 格 非常 细密 时 ， 离 散 方程 才 满 足 积分 守恒 ; 

图 5. 1 为 计算 控制 体 图 。 大 写字 母 为 控制 体 的 中 心 点 ，E、 由 表示 东西 方向 上 控制 
体 中 心 点 , N、S 表示 南北 方向 上 控制 体 中心 点 ;U、D 表示 上 下 方向 上 控制 体 中 心 点 。 e, 
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v 表示 东西 方向 界面 ，n，s 表示 南北 方向 界面 ，u，d 表示 上 下 方向 界面 。 


UW | U UE 
^4 
US 
NE 


‘ 






DW 


图 5. 1 计算 控制 体 图 
在 以 P 点 为 中 心 的 有 限 体积 上 对 控制 方程 进行 积分 ， 经 过 宛 繁 的 推导 ， 最 后 得 出 全 
隐 式 时 间 积 分 方案 下 的 三 维 瞬 态 对 流 - 扩 散 问题 的 离散 方程 为 : 


Andy = Ay Gy + Arby +Asps +A by o Au), + Ardy +S (5. 1) 
离散 方程 的 通用 表达 式 
Arby = NU En F] g (5.2) 
nb-E At 


式 中 ， 下 标 nb 表示 相 邻 节点 。 对 于 三 维 问题 ， 相 邻 节点 包括 W、E、S、N、U 和 D。 
上 标 “0” 表 示 上 一 时 层 。 
APS, 为 源 项 ， 其 表达 式 为 


S, = JSAV + Avebye + Ay Pw + AseGse + Ayghug 
+ Ayr dys + Auw Puy + ÁpsÓpg + Apw Pow + Aun duy (5. 3) 
+ Aysdys + Ápy doy + AnsPns 


Ap» Ays Aug FEIRA: 


P [p 
f nb=E Zi P 


Ag = D,AQP, |) + []-5,.0]] Ay = Dy AQP. ow) + [| F,,0] 
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Ay = D,AGP,,|) + [I-5.,. 0] , As = D,A(P.,. * [5.9] 


A, = Dy AdE D [[-5..4l] » A, = DAQPLD + [I-&..9]] 
Ay = Dan *Du Ay = -Dan - Dir, 
Ag =~Dy,- Ds Asy "D D, 
Ag "D, *tD Aw = -Dy -Dav 
Ay = -Daa -Diet Ay = Dj, + D3, 
Ay "Do, + Dy, Ag = Ds, - Drs, 
Ay = -Doa Dy,» Aps = Di, + D, 
| D, = Dy, + Dy, - Dy, + Dy -D,,» D, = Diw ~ Dan +D, -Doi + Dry 
D, = Dy, + Dy, -Diaw + Dy, Do,» D, = Dy, ~ Dy, + Dy, -Dou + Do, 
D, = Dozu + Dise -Diaw + Dz, ~Da» D, D, -Ds + Diaw -Dan + Das, 
F, = pU*AS,, F, = pU?AS,, F, = pU°AS, 


(436) aiAS 不对，j 求 和 ，6, 为 克 罗 尼 柯 尔 记号 ) 


一 一 一 一 -一 一 


”4 





P, (Peclet 数 ) 


以 上 式 子 ，F 表示 通过 界面 上 单位 面积 的 对 流质 量 流量 ，D 表示 界面 的 扩散 传导 性 


不 同 格式 的 区 别 仅 在 于 AP) KHARE, A(R ERUR 5. 2 Pra 
表 5.2 差分 格式 
alz.) 
o -osp 


l FP, 
指数 格式 HETE 


乘 方 格式 il. (1-0:118]) 


#52 
符号 “| | ”表示 取 各 量 中 之 最 大 值 。 
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5. 4 小 结 


本 章 围绕 如 何 数 信 求 解 控制 方程 ， 对 比 了 离散 方法 ， 讨 论 了 离散 方法 的 选取 和 本 文 
采用 的 离散 策略 、 物 理 量 的 配置 方法 。 重 点 讨论 了 如 何 计算 控制 体 及 其 控制 面 的 几何 要 
素 ， 如 何在 直角 坐标 系 下 离散 控制 方程 。 
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第 六 章 ”离散 化 方程 的 求解 算法 


6 1 ” 非 线性 方程 的 线性 化 


经 过 第 五 章 方程 的 数值 离散 处 理 ， 在 决定 方程 的 哪些 项 采用 隐 式 、 哪 些 项 采用 显 式 
后 ， 形 成 了 方程 的 系数 和 源 项 ， 将 非 线性 的 RANS 微分 方程 转化 为 形式 上 为 线性 的 代数 
方程 。 对 于 每 个 控制 体 ， 动 量 方程 、k 方 程 、 方程 都 可 以 写成 统一 的 离散 形式 swm; 

Arby = Ay ty +Arbe + Aghs + Aupy + Anha + Arty +S (6.1) 


APP 是 待 求 物理 量 ，4 表示 系数 ，$ 表示 源 项 ， 下 标 表示 节点 编号 ， 如 果 把 所 有 
控制 体 的 守恒 离 散 方程 联 立 起 来 ， 得 到 一 个 线性 方程 组 *****”: 


OR RR (6.2) 


A, p kf ka + Á, p kadr ra 


6 2 ”求解 代数 方程 组 的 两 种 方法 中 >** 各 名 


6. 2. 1 直接 解法 


所 谓 直接 解法 (direct method) 是 指 通过 有 限 步 的 数值 计算 可 以 获得 代数 方程 真 解 . 
的 方法 。 其 中 ， 高 斯 消 元 法 是 一 类 用 得 较 多 、 具 有 普遍 意义 的 方法 。 用 高 斯 消 元 法 求解 
多 维 问题 的 代数 方程 时 ， 先 要 将 方程 的 系数 通过 消 元 化 为 上 三 角 和 矩阵 , 然后 从 上 三 角形 
的 下 方 顶 点 开始 ， 逐 一 回 代 计算 出 各 个 未 知 量 。 由 于 高 斯 消 元 法 系数 的 存储 是 方程 个 数 
的 平方 ， 同 时 所 用 乘法 的 次 数 近似 正比 于 方程 个 数 的 三 次 方 ， 因 此 当 方 程 个 数 较 大 时 ， 
计算 所 需 的 时 间 和 内 存 都 会 急剧 增加 。 有 人 做 过 试验 :对 于 二 维 正 方形 区 域内 的 拉 普 拉 
斯 方程 ， 当 采用 3200 个 节点 时 ， 采 用 高 斯 消 元 法 所 需 的 计算 次 数 是 完成 一 次 雅 可 比 迁 
代 所 需 计算 次 数 的 168 倍 。 所 以 当 求解 的 节点 数 较 多 时 ， 高 斯 消 元 法 的 效率 非常 低 ， 因 
而 限制 了 这 种 方法 在 粘性 方程 求解 中 的 应 用 。 
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6. 2. 2 RM 


对 于 本 文 所 讨论 的 问题 ， 从 形式 上 而 言 ， 尽 管 所 需求 解 的 是 线性 代数 方程 ， 但 是 其 
实质 是 求解 非 线 性 的 NS 方程 。 离 散 方程 的 系数 、 源 项 都 是 用 上 一 个 选 代 时 间 步 的 物理 
量 来 表达 的 ， 都 是 最 终 待 求 未 知 量 的 函数 。 这 样 ， 整 个 问题 的 求解 必然 具有 是 迭代 性 质 
的 : 即 先 假定 一 个 初始 的 状态 场 (具体 编程 实现 是 计算 起 始 时 ， 需 要 调用 初始 化 函数 ) ， 
根据 此 计算 出 离散 方程 的 系数 ， 然 后 求解 方程 得 到 改进 值 ， 进 而 重新 计算 方程 的 系数 ， 
再 求解 ，…… ， 如 此 反复 ， 直 到 方程 的 求解 达到 收敛 。 

由 此 可 见 ， 方 程 的 系数 都 是 临时 的 ， 如 果 采 用 直接 求解 ， 得 到 精确 值 也 只 是 这 样 一 
组 临时 系数 的 解 而 己 。 既 然 方程 的 系数 本 身 就 是 有 待 改进 的 ， 把 相应 的 精确 解 计算 出 来 
既是 费时 、 费 内 存 的 ， 也 是 没有 必要 的 . 采用 从 代 法 计算 时 ， 通 过 添加 判 据 减 小 或 终止 
方程 迭代 ， 以 在 改进 代数 方程 的 系数 后 再 进一步 求解 。 因 而 对 于 非 线 性 问题 的 线 件 化 方 
程 的 求解 ， 采 用 迭代 法 是 很 有 必要 的 ， 也 是 经 济 的 。 

当前 ， 在 数值 计算 领域 ， 用 于 求解 线性 方程 组 的 迭代 方法 很 多 。 比 较 传统 的 有 雅 可 
比 和 迭代 法 、 高 斯 一 赛 德 迭 代 法 。 雅 可 比 和 迭代 法 是 把 任何 一 点 未 知 数 的 更 新 使 用 上 一 轮 和 迭 
代 中 所 获得 的 各 相 邻 点 的 值 来 计算 ， 其 迭代 前 进 方向 并 不 影响 迭代 收敛 的 速度 ;高 斯 赛 
德 运 代 法 是 每 一 步 的 计算 总 是 取 相 邻 点 的 最 新 值 来 进行 ,边界 条 件 的 影响 传 入 区 域内 部 
的 快慢 是 不 同 的 ， 和 迭代 前 进 方向 对 收敛 性 有 较 大 影响 。 在 计算 流体 力学 领域 ， 比 较 受 
研究 者 欢迎 的 迭代 法 有 超 松弛 迭代 法 〈 亚 松弛 迭代 法 ) 、 交 蔡 隐 式 和 迭代 法 、 强 隐 式 迭代 


6. 3 ” 超 松 弛 迭代 法 (SOR) 、 亚 松弛 迭代 法 (SUR) 


该 方法 是 把 任何 一 点 未 知 数 的 更 新 p*“ 使 用 上 一 轮 选 代 中 所 获得 的 各 相 邻 点 的 什 

J 和 本 轮 选 代 中 所 获得 的 最 新 值 g” 加 权 平 均 来 计算 中 中 : 
pg alg” -9") (6. 3) 
At, a 被 称 为 松弛 系数 ，0<a s2 。 其 中 ，a -0 时 ， 是 雅 可 比 迭 代 ，c =1 时 ， 
RN PS: 0<a <1 时 ， 是 亚 松弛 园 代 (CSUR) ; 1<a <2 时 ， 是 超 松弛 迭代 
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(SOR) 。 
从 式 (5.4) 可 知 ， 当 相 邻 两 轮 的 迭代 值 永远 具有 相同 的 正 负 号 时 ， 采 用 超 松弛 迭 代 
可 以 加 快 收敛 速度 。 对 于 非 线 性 问题 所 形成 的 代数 方程 ， 多 采用 亚 松弛 和 迭代 法 ， 以 降低 
待 求 量 的 变化 率 ， 避 免 迁 代 发 散 。 


6.4 simple 算法 


求解 离散 后 的 动量 方程 时 ， 由 于 压力 场 和 速度 场 均 是 未 知 量 , 而且 没有 关于 压力 本 
身 的 控制 方程 ， 压 力 和 速度 的 关系 只 是 隐 含 在 连续 性 方程 中 。 

为 此 Patanker 和 Spalding 提出 Simple 算法 (Semi-Implicit Method for 
Pressure Linked Equation)， 其 基本 思想 是 ““ 趾 ， 如 果 压 力 场 是 正确 的 那么 根据 该 压 
力 场 求解 动量 方程 得 到 的 速度 场 必定 是 正确 的 满足 连续 性 方程 ， 为 此 ， 可 在 压力 场 未 
知 的 情况 下 ， 先 假定 一 个 压力 场 ， 该 压力 场 一 般 不 可 能 就 是 所 求 之 值 ， 因 而 它 得 到 的 速 
度 场 也 不 会 正确 ， 不 会 满足 连续 性 方程 ， 但 是 可 用 连续 性 方程 对 压力 进行 修正 ， 得 到 压 
力 修 正方 程 ， 进 而 对 速度 进行 修正 ， 进 一 步 求解 能 满足 连续 性 方程 的 速度 场 ， 然 后 以 修 
正 后 的 速度 场 计算 动量 离散 方程 的 各 个 系数 ,修正 后 ,在 压力 场 的 基础 上 进行 新 一 次 的 
计算 ， 如 此 ， 直 到 代入 动量 方程 的 压力 场 求 出 的 速度 场 能 够 得 到 连续 性 方程 为 止 ， 因 此 
问题 的 关键 是 得 到 压力 修正 方程 。 

压力 修正 法 是 目前 工程 上 使 用 最 为 广泛 的 流 场 数 值 计算 方法 ， 其 实质 是 迭代 法 。 在 
每 一 时 间 步 长 的 运算 中 ， 先 给 出 压力 场 的 初始 猜测 值 ， 据 此 求 出 猜测 的 速度 场 ， 再 求解 
根据 连续 方程 导出 的 压力 修正 方程 ， 对 猜测 的 压力 场 和 速度 场 进行 修正 。 

BADR: 

步骤 (1) : 给 定 初始 速度 场 ， 对 可 压缩 流动 还 要 给 定 初始 压力 场 和 温度 场 

步骤 (2) ;求解 动量 方程 得 到 本 次 迭代 的 速度 uw，v，w 

步骤 (3) ; 求解 压力 修正 方程 得 到 压力 修正 量 p, Hp 修正 压力 场 和 速度 修 
IER wu, y , w ， 并 且 改 进 速度 场 。 

步骤 (4) : 求解 上 方程 和 e 方程 ， 进 而 得 到 jp 和 离散 方程 的 系数 。 

步骤 (5) ， 对 可 压缩 流动 求解 能 量 方程 得 到 温度 场 
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步骤 (6) : 如 果 达 到 收敛 指标 则 停止 计算 , 否则 返回 到 步骤 (2) 进行 下 一 次 迭代 计算 。 


6.4.1 SIMPLE 算法 基本 原理 “所 入“ 名 


流 过 每 个 单元 面 上 的 对 流通 量 根据 “猜测 ”速度 估算 。 首 先 使 用 一 个 猜测 的 压力 场 
来 解 动量 方程 ， 得 到 速度 场 ， 接 着 求解 通过 连续 方程 建立 的 压力 修正 方程 ， 得 到 压力 场 
的 修正 值 ， 然 后 利用 压力 修正 值 更 新 速度 场 和 压力 场 ; 检查 结 果 是 否 收 伐 ， 若 不 收敛 ， 
以 得 到 的 压力 场 作为 新 的 猜测 压力 场 ， 重 复 该 过 程 。 

为 了 启动 该 迭代 过 程 ， 需 要 提供 初始 的 ， 带 有 猜 测 性 的 压力 场 与 速度 场 。 随 着 迭代 
的 进行 , 这 些 猜测 的 压力 场 与 速度 场 不 断 改 善 , 所 得 到 的 压力 与 速度 分 量 逐 渐 逼 近 真 解 。 
假定 压力 的 应 力 初始 猜测 值 p*, AF I, J, K 代表 控制 体积 节点 编号 , i, j,k 
代表 控制 体积 界面 编号 ， 则 有 




















Payk 7 pis ' uma + Pray a) (6. 4) 
Box = iss = FA + Park ies | (6. 5) 
Pije = los : 7 o. + Prag A (6. 6) 
Prax = Pisa + Dryax a (6. 7) 
ee piss Aiet s piyar ica (6.8) 
Pisin 7 区 Do Prax a (6. 9) 


动量 方程 如 下 : 
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s div(pVu) » div(u gradu)-2 +5, 
x 





Ee aure) -nugredv) -P+s, (6. 10) 





(ow) + div(pVw) = div( grad w)-2 45, 
ot oz 
在 非 结构 网 格 上 的 动量 离散 方程 恕 下 : 
x 
6, y xHy p k ™ p GS, + (Px -pix )AvAz +b 
Arsis” Yaara +( Pix - Pi jux )Axz +b (6. 11) 
d, P x Wi x ™ Saswat +(Prss - pua) AxAy +b 


LP, HE e RATA PISA Par Puo Pus Pas Pujas? Pissa 
代表 界面 上 e 的 压力 ， 系 数 Ar ，Ay ，Az 是 界面 e 的 终点 与 起 点 的 坐标 之 差 ， 它 是 有 正 
负 的 。 各 界面 的 起 点 及 终点 次 序 按 逆 时 针 排列 。 

一 .速度 修正 方程 的 建立 | 
根据 方程 组 (6-10) 得 出 网 格 上 的 界面 速度 方程 (及 动量 插值 方程 ) 为 ， 








ae [att |. Bava Fae) (6. 12) 
uy M yg E Pw | ; (Pox -Fssx) (6. 13) 
Viu 5 pue (e) (Baas -B x) (6. 14) 
Vix zs 山寺 7 cae R (Box -B ax) (6. 15) 





n A, rk 
"atus n (Bisan Bk) (6. 16) 
I LJ, k 
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~ A 

1J,K 

Wisaa W rka jy (Bx -Bra) 
13K Ji, yi 





(6. 17) 


其 中 x ，v ，w 称 为 伪 速 度 ， 界 面 上 的 物理 量 w ，y ，w 通过 控制 体积 P 与 控 
制 体 积 E 中 心 节点 处 的 物理 量 线性 插值 得 出 。 





N Ny ; Xx 
= 2s uk Die Sk - 2j uk 4 
u 2 zl £A NE Jia a | Ji 7h (6 18) 
WK z b : 
Gk yg Xu Gk Ma 
IK IJK TK 
X % n 
Dito +b yx whe thy x Die +b yk 
x54 .*7X 
TE : A 2% 6.19) 
Gig asx 44a "ak Ñ -Xa 
iJK IK TAJK 
N, N, 
. PT +b, ,x goin Sata *b, ak yey 
Vine = || =E os | Fa uL 2i. (6. 20) 
Gik Ji Yra Q, 7-LK yr Yin 
EJK ISAK 
N. Ns 
anu, +D, AU, +5, 5, 
s 2: > DN PE a Viet + 2 ie iua ja 7 Y; (6. 21) 
Lp * 
uA. Yau 7 Ys BAK Yiu ~ Y; 
IJK TAK 
N, N, 
M Data +b, 5% zas Sasis 4855-4 2 
wisa = || = KE C. (6. 22) 
d, rk 三 一 zx- d, pk Zr 一 Zr 
IJK iJ Ki 
l N N, 
* 
" D Culla tb, yx gh Data +h kn dius 
Wi ka = || E Kal Westy] ESO 00 Ck Tk (6. 23) 
d, jk Zga 77x 8r ka Zkn 77k 
DLE iJ K+ 
《 


式 (6.18) — (6.23) 中 aj jk ，biy x 的 表达 式 如 下 : 
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N AV 
Ayk" 2^ + Por -$, AV 
= : (6. 24) 
u,, AV 
bjx _ (pmsxAr) sae *$,AV 


RR FE NE BIE p* ， 即 (6. 4) — (6.9) 所 得 到 的 p. (6.12) — 
(6.17) 可 以 求 出 网 格 边界 上 的 猜想 速度 。 

然后 进行 速度 修正 方程 的 建立 , 假定 压力 修正 值 P' 通 过 压力 修正 方程 (6-39) 已 经 
求 出 ， 引 入 SIMPLE 算法 中 略 去 邻 点 速度 修正 值 的 思想 ， 得 到 界面 上 的 速度 修正 方程 如 
下 : 

















A, 
unk zt [Bos z Bix) (6. 25) 
Gx jx 
j 4 ; i 
Wisk -| < (Pox -Px) (6. 26) 
rk iar 
D A j " " 
Vij = - [PR = P ik ) (6. 27) 
"n hjx 
t A 1 
Vi, jak -| tie | (Fx -Pax (6. 28) 
RR ME 
f A ; : 
Was” a (Bs -Px ) (6. 29) 
9p pk IJk 
; A i : 
Waa” dd [Es -B xa] (6. 30) 
KEZ IJke 


然后 ， 把 (6.12) — (6.17) PGB ues lase’ Vx Vra Work? 


VETE Vo. wy w, 加 上 (6.26) — (6.30) Br 


mJ ? Vx Leja IJA i Idka 
FEA ue Mare? V Vipax? Woe? Was 最 终 得 到 控制 体积 界面 上 的 “ 正 
确 ” 速 度 方程 : 
即 U-U'«U (6.31) 
iX HU Riu, v, w 
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在 这 里 分 别 表 示 为 
A 
Wx tl) (Fase Pe) (6.32) 
EN. id, K 
P A 
las “Uy cm [Ere Bae] (6.33) 
1J,K iJ, K 
e A ; : 
wire | E (Boax- Box) (6.34) 
rk jinx 
. fA 
Vra TV, jr end (Bx TÉLIAK (6. 35) 
Lje K 
« A , ' 
Wisp mW, e| oe (Bria Bir) (6. 36) 
22 P 
* A + , 
Wi pa mW, en (Bx - Box) (6. 37) 
IJ. k 


控制 体积 节点 上 正确 的 速度 场 同 理 也 可 以 求 出 ， 在 这 里 不 作 推 导 。 
Z. 压力 修正 方程 的 建立 
散 度 形式 的 连续 方程 如 下 : 


ap , (ou) (pv) alow) o (6. 38) 
of a ay az 


在 时 间 间 隔 At 内 对 控制 体积 P 作 积分 ， Bp EP pyt, 采用 全 隐 式 格式 ， 


可 得 : 





BaP AV + [(puA), - (pua), ] + [(oud), - (pua), ] + [(pu4), -(pud),]=0 (6.39) 


将 正确 的 速度 值 ， 即 〈6. 32) -— (6.37) 所 得 到 的 wj， laro Vr Vra 
mw jw yi 代入 连续 离散 方程 (6.39) 中 ， 得 到 压力 修正 方程 : 
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Qa KPig kK = Fas x Prag tay Pis. + Ug KP rsa 


+r | axPrsax ty xubrs ga * T1 kar) aa D (6. 40) 
式 中 ， 系 数 

Graze = (PdA) jk (6. 41) 

Gas = (pd4), ,x (6. 42) 

a, jx = (pdA), jx (6. 43) 

2, jx = (pdA), x (6. 44) 

4, , xn" (0d4), | i (6.45) 

8, jk. = (pd^), , , l (6.46) 

(6. 47) 


Gisk "Opa pg lras tarsar tU aK tA KA ts 


bsx m (ouA) e -(ou'A) y +(pv'A) y -(pv'A), x +(owA) - (ow) pa 
à Prox ~ Pria AV 
At 
(6. 48) 
从 方程 (6. 39》 可 解 出 节点 上 的 压力 修正 值 P', x ， 再 通过 线性 插值 可 得 出 界面 上 
的 压力 修正 值 P。 按 下 式 得 到 修正 后 的 压力 场 : 


P Pu (6. 49) 


6.4.2 ” 选 代 前 的 亚 松 弛 预 处 理 


本 文采 用 simple 迭代 算法 进行 方程 求解 时 ， 压 强 修正 方程 和 速度 修正 方程 容易 发 
散 ， 为 了 提高 求解 的 稳定 性 ， 除 非 对 方程 进行 低 松 弛 处 理 。 同 时 ， 为 了 节约 内 存 存放 所 
有 控制 体 物理 量 的 中 间 结 果 ， 对 式 〈6. 31》 和 (6. 49) 统一 进行 如 下 处 理 外 ; 
9" =a" «(1- a)g" (6.52) 
XP a 为 低 松弛 因子 ， 修 正 量 大 到 足以 不 断 修正 ， 但 又 应 该 足够 小 ， 以 保证 稳定 计 
算 ， 所 以 低 松弛 因子 x 一 般 取 在 1 和 0 之 间 。 
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O° 表示 方程 组 的 迭代 满足 收敛 条 件 的 求 出 值 ， 是 方程 的 真实 根 ， 它 是 满足 每 个 
控制 方程 的 。 将 加" 代替 式 (6. 1) 中 的 g 得 到 : 


5 
có, = s +5, (6. 53) 


式 中 ，i = 0,1...5 表示 周围 的 六 个 相 邻 节点 。 
将 式 (6. 52) 代入 式 (6. 53) 左边 ， 碳 边 结 合 高 斯 - 赛 德 夫 代 法 的 思想 ， 在 处 理 相 邻 点 
时 使 用 最 新 值 ， 得 到 N+1 步 选 代 关 系 式 : 


C, NA 1- x +1 
“gs -Cg + Seg +s, (6. 54) 
从 根本 上 说 ， 正 确 选 择 低 松弛 因子 ， 会 有 效 的 降低 了 成 本 
在 计算 时 ， 低 松弛 因子 对 动量 方程 、 连 续 性 方程 、xk JEA e 方程 的 取 值 不 同 ， 如 
下 表 所 示 : 


表 6.1; 松弛 因子 


CE CL e 







关于 方程 的 最 大 和 迭代 步 数 ， 在 迭代 开始 ， 取 值 小 一 些 ， 增 强 方程 的 稳定 性 ， 在 求解 
后 期 ， 取 值 大 一 些 ( 细 算 ) ， 提 高 方程 的 求解 的 效率 。 本 文 按照 表 5. 2 取 值 : 
R 6. 2: 方 程 最 大 迭代 步 


动量 方程 连续 方程 Kk 方程 2 方程 


选 代 的 收敛 标准 精度 本 文 取 : 小 于 1E-4， 

对 于 选 代 的 收敛 精度 ， 本 文 取 相对 误差 ， 分 别 以 各 个 代数 方程 的 相对 残 差 为 准则 : 
a 

R, =-= 


f Norm, 
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Ab, RE > 为 迭代 完成 时 ， 每 个 控制 体 方程 的 左右 两 边 的 差 。Normz 为 相对 误差 
参考 值 。 


6. 6 ”求解 器 设计 


求解 过 程 从 几何 要 素 初始 化 、 流 场 初始 化 开始 ， 本 文 取 入 流 边界 值 作为 全 场 压 力 ， 
速度 以 及 t 和 的 初始 值 ， 并 且 强 制 取 壁面 边界 上 的 速度 为 零 。 求 解 中 需要 夫 代 求解 线 
性 方程 组 ， 本 文 先 计算 出 动量 方程 的 系数 ， 然 后 按照 山  W 的 顺序 ， 依 次 求解 出 各 
方向 的 速度 场 后 ， 再 作 压力 场 推进 。 最 后 依次 计算 上 和 输送 方程 的 系数 ， 先 后 对 这 两 
个 方程 进行 求解 。 具 体 过 程 见 图 6.1 

结合 图 6. 1， 本 文采 用 标准 FORTRAN95 语言 ， 设 计 了 RANS 方程 求解 器 程序 。 


6. 7 本章 小 结 
本 章 就 RANS 方程 线性 化 离散 后 得 到 的 代数 方程 组 的 求解 进行 了 讨论 ， 选 用 了 
SIMPLE 算法 求解 方程 , 并 对 三 维 SIMPLE 算法 进行 了 详细 介绍 。 为 了 加 强求 解 的 稳定 性 ， 


在 方程 每 步 和 迭代 之 前 ， 对 线性 方程 进行 了 亚 松弛 预 处 理 。 然 后 就 方程 组 求解 的 收敛 性 进 
行 了 判断 ， 给 出 了 收敛 性 条 件 。 最 后 结合 具体 求解 过 程 的 分 析 ， 设 计 了 方程 求解 器 。 
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开始 


读 入 网 格 数据 ， 凡 何 要 素 初始 化 


读 入 控制 文件 流 场 初始 化 


计算 动量 方程 和 各 线性 方程 的 系数 ， 并 依次 求解 fd 


U, V. WHE 
e | es 
求解 器 





连续 性 方程 压力 场 计算 [M 
P| mn 
计算 x 方程 的 系数 ， 并 求解 方程 E 新 节点 值 
计算 e 方程 的 系数 ， 并 求解 方程 EMESA 
FEB EINE 
代 收 敛 条 件 
= 
是 
输出 结果 数据 ， 数 据 文件 ， 绘 图 文 
件 


6.1 求解 器 设计 流程 图 
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第 七 章 算 例 及 结果 分 析 


为 了 验证 方程 求解 器 的 稳定 性 、 正 确 性 , 同时 方便 数据 的 检查 和 计算 程序 调试 ， 本 
文 首先 进行 了 亚 音速 流 计算 ， 然 后 对 超 音 速 流 计算 ， 最 后 计算 结果 分 别 和 风 洞 数据 进行 
对 比 。 些 时 ， 求 解 是 基于 高 速 飞行 器 的 ， 在 固 壁 边界 上 仅仅 满足 无 滑 移 条 件 即 可 ， 大 大 
降低 了 程序 处 理 的 复杂 性 。 


711 Mk SHR" 


细 长 旋 成 导弹 体 ,长 细 比 为 10, 圆柱 长 为 80, 弹头 长 为 20, 口 径 为 10, 在 攻 角 a =0， 
在 1000x150 的 流 场 中 的 速度 , 温度 , 压力 分 布 情况 , 这 里 仅 反 映 对 称 面 流 场 。 这 里 边界 条 
FE: t 界面 为 入 流 边 界 条 件 U,，U, -U, =0，-x 界面 为 出 流 边界 条 件 ， 导弹 体 周 
边 流 场 定义 为 远 场 条 件 ， 对 称 面 定 义 为 对 称 边 界 条 件 

来 流速 度 为 (为 0.6 马赫 ， 在 上 面 的 条 件 下 所 得 到 的 残 差 收敛 史 , 压力 分 布 图 ， 
速度 分 布 图 ， 温 度 分 布 图 ， 速 度 和 撩 量 图 ， 压 力 分 布 曲 线 图 ， 速 度 分 布 有 曲线 图 ， 温 度 分 布 
曲线 图 ， 分 别 如 图 7. 1，7.2，7.3，7. 4，7. 5, 7. 6, 7.7 所 示 (只 显示 对 称 面部 分 ): 


Residuals 
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图 7.1 细 长 旋 成 导弹 体 对 称 面 上 残 差 收敛 史 
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ECG 


2 细 长 旋 成 导弹 体 对 称 面 上 的 压力 分 布 云图 








图 7.3 细 长 旋 成 导弹 体 对 称 面 上 的 马赫 数 分 布 云图 
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图 7.4 细 长 旋 成 导弹 体 对 称 面 上 的 温度 分 布 云图 
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sx 





图 7.5 细 长 旋 成 导弹 体 尾部 速度 矢量 图 
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图 7.6 细 长 旋 成 导弹 体 对 称 面 速度 分 布 曲 线 


Static 
Pressure 
(pascal) 





Position (mm) 


图 7.7 细 长 旋 成 导弹 体 对 称 面 压力 分 布 曲 线 
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图 7.8 细 长 旋 成 导弹 体 对 称 面 温度 分 布 曲 线 


来 流速 度 为 0 为 2 马赫 在 上 面 的 条 件 下 所 得 到 的 残 差 收敛 史 , 压力 分 布 图 ， 速 度 分 
DA, 温度 分 布 图 , 速度 矢量 图 ， 压 力 分 布 曲线 图 ,速度 分 布 曲线 图 ,温度 分 布 曲线 图 ， 
分 别 如 图 7. 9, 7. 10, 7. 11, 7.12, 7. 13, 7. 14, 7. 15, 7. 16 所 示 《〈 只 显示 对 称 面部 分 )， 





tterations 


图 7.9 细 长 旋 成 导弹 体 对 称 面 上 残 差 收敛 史 
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图 7.10 细 长 旋 成 导弹 体 对 称 面 讨 力 分 布 云图 
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图 7. 11 SCRE SE HORE SIR DS BC a ER 
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图 7.12 细 长 旋 成 导弹 体 对 称 面 温度 分 布 云图 
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图 7.13 细 长 旋 成 导弹 体 尾部 速度 矢量 图 
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图 7.14 细 长 旋 成 导弹 体 壁 面 速 度 分 布 曲 线 


Static 
Pressure 
(pascal) 





Position (mm) 


Bd 7.15 细 长 旋 成 导弹 体 壁 面 压力 分 布 曲线 
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图 7.16 绍 长 旋 成 导弹 体 壁面 温度 分 布 曲线 


来 流速 度 为 U0: 为 不 同 马 赫 数 时 , 阻力 系数 计算 值 与 实验 值 如 表 7. 1 所 示 ， 阻 力 系 数 
计算 值 和 实验 值 比较 如 图 7. 17 所 示 | 


表 7.1 阻力 系数 计算 值 与 实验 值 及 误差 


阻力 系数 计算 值 | 阻力 系数 实验 值 ”| 相对 误差 











«m 
[ ss | usa l uuB ux 
[sw | wm | wmm 8 
o a j o 
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0. 3012 0.3119 3. 4% 
—o—it 8 ff 
一 一 实验 值 





n EA 
图 7.17 细 长 旋 成 导弹 体 阻力 系数 变化 对 比 曲线 
从 图 7.17 可 以 看 出 不 管 在 亚 音速 流 时 还 是 在 超 音速 流 时 ， 对 细 长 旋 成 导弹 体 计 算 
的 结果 基本 还 是 准确 的 。 阻 力 系数 都 是 在 1. 5 马赫 时 达到 最 大 。 


7. 2 Ree 


未 制导 弹 体 ， 长 细 比 为 8. 5， 圆 柱 长 为 80, 口径 为 10， 强 头 为 圆 形 ， 弹 头 半 径 为 5， 
在 攻 角 a =0， 在 1000x150 的 流 场 中 的 速度 , 温度 , 压力 分 布 情况 , 这 里 仅 反映 对 称 面 流 
场 。 这 里 边界 条 件 是 : tx 界面 为 入 流 边界 条 件 U,，U, =U, =0，-x 边界 为 出 流 边界 条 
件 ， 导弹 体 周边 流 场 定义 为 远 场 条 件 ， 对 称 面 定义 为 对 称 边界 条 件 

来 流速 度 为 0: 为 0.6 马赫， 在 上 面 的 条 件 下 所 得 到 的 残 差 收敛 史 , 压力 分 布 图 ， 速 
度 分 布 图 温度 分 布 图 ， 压 力 分 布 曲 线 图 ， 速 度 分 布 曲 线 图 ， 温度 分 布 曲 线 图 ， WER 
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量 图 , 分 别 如 图 7. 18, 7. 19, 7. 20, 7. 21, 7. 22, 7. 23, 7. 24, 7. 25 所 示 ( 只 显示 对 称 面部 分 ): 
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7.18 未 制导 弹 体 对 称 面 上 残 差 收敛 史 
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图 7. 19. 末 制 导 阐 体 对 称 面 上 速度 分 布 云图 
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图 7. 20. 末 制 导弹 体 对 称 面 上 压力 分 布 云图 
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图 7. 21. 未 制导 弹 体 对 称 面 上 温度 分 布 云 图 
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图 7.22 末 制 导弹 体 对 称 面 上 速度 曲线 分 布 





2.000005 
1.004605 
0.004000 
-200 9 200 400 600 $00 4003124608 
Position (mm) 


图 7. 23 末 制 导弹 体 对 称 面 上 压力 分 布 曲 线 图 
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图 7. 24 末 制 导弹 体 对 称 面 上 温度 分 布 曲 线 图 
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图 7. 25 “未 制导 弹 体 对 称 面 上 尾部 速度 矢量 分 布 图 

来 流速 度 为 中 为 2.0 马赫 ， 在 上 面 的 条 件 下 所 得 到 的 残 差 收敛 史 , 压力 分 布 图 ， 速 

度 分 布 图 ， 温 度 分 布 图 ， 压 力 分 布 曲 线 图 ， 速 度 分 布 曲线 图 ， 温 度 分 布 曲线 图 , 速度 和 

BA, 分 别 如 图 7. 26, 图 7.27， 图 7.28， 图 7.29， 图 7.30， 图 7.31， 图 7.32， 图 7.33 
BUR (只 显示 对 称 面部 分 ); 
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图 7. 26 末 制 导弹 体 对 称 面 上 残 差 收敛 史 


图 7. 27 末 制 导弹 体 对 称 面 上 速度 分 布 云图 
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图 7. 28 末 制 导弹 体 对 称 面 上 压力 分 布 云图 





”图 7.29 末 制 导弹 体 对 称 面 上 温度 分 布 云图 
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7. 30. 末 制 导弹 体 对 称 面 上 速度 分 布 曲线 图 
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图 7.31 末 制 导弹 体 对 称 面 上 压力 分 布 曲线 图 
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图 7. 33. 末 制 导弹 体 对 称 面 上 尾部 速度 矢量 图 
来 流速 度 为 0: 为 不 同 马琳 数 时 ， 阻 力 系数 计算 值 和 实验 值 如 表 7.2 BUR, MHA 
数 变 化 曲线 图 如 图 7. 34 所 示 
表 7.2 阻力 系数 计算 值 和 实验 值 及 误差 
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7.34 末 制 导弹 体 阻力 系数 变化 对 比 曲 线 


71. 3 GRRE 


从 表 7. 1 图 7.17 和 表 7.2, 图 7.34 可 以 看 出 ,阻力 系数 基本 吻合 , 从 图 7.13 和 7.3 
可 以 看 出 , 在 两 种 类 型 导弹 体 尾 部 都 出 现 涡流 , 从 图 7. 28， 图 7. 29， 图 7. 30 可 以 看 出 在 
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导弹 体 弹头 部 分 有 明显 的 激 波 现象 , 而 从 图 7.2， 图 7.3， 图 7. 4 在 导弹 体 弹 头 部 分 却 看 
不 出 有 上 明显 的 激 波 现象 , 这 些 图 象 、 曲 线 都 正确 反映 了 高 速 飞行 器 流 场 的 一 些 基 本 特征 ， 
对 高 速 飞行 器 空气 动力 学 数值 分 析 有 一 定 的 意义 。 


7. 4 小结 


为 了 检验 求解 器 的 准确 性 , 本 章 对 超 音速 和 亚 音速 流 分 别 进 行 数值 计算 , 并 给 出 求 
解 结果 , 计算 结果 表明 ， 不 管 在 亚 音速 流 计算 还 是 在 超 音速 流 计算 者 是 比较 准确 的 。 


59 


中 北大 学 学 位 论文 


结论 


本 文 讨论 了 三 维 洪流 流动 的 计算 理论 ， 从 三 维 直角 坐标 系 下 的 流动 控制 方程 出 发 ， 
介绍 了 直角 坐标 系 下 的 三 维 流动 控制 方程 ， 并 给 出 了 计算 控制 体 中 的 离散 方程 。 

对 各 种 洪流 模型 进行 详细 对 比 ， 采 用 上 ~ 满 流 模 型 进行 潍 流 的 数值 模拟 。 

详细 介绍 了 三 维 SIMPLE 算法 ， 并 应 用 三 维 SIMPLE 算法 进行 不 可 压 流 场 计算 。 编制 
了 三 维 淇 流 流 场 计算 的 相应 软件 , 并 对 细 长 旋 成 导弹 体 和 未 制导 弹 体 的 流动 进行 了 数值 
HH, 成功 地 模拟 出 流动 现象 。 

由 于 编制 和 调试 三 维 满 流 流 场 计算 程序 工作 量 很 大 ， 和 需要 的 知识 很 多 . 尤其 对 每 一 
个 算 例 的 调试 是 很 费时 的 ， 所 以 还 需要 进行 更 多 三 维 流 场 数值 模拟 ， 以 进一步 检验 和 修 
正本 文 的 算法 。 
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学 的 发 展 而 发 展 的 . 从 较 早 的 从 宏观 角度 出 发 求解 Euler 方 程 , N-S 27 8, 到 考虑 真实 气体 影响 带 化 学 反应 和 电离 的 N-S 方 程 的 求解, 以 及 从 微观 角度 出 发 用 蒙特 

卡 洛 法 对 气体 分 子 运 动 直 接 模拟 求解, 即 代表 计算 动力 学 理论 的 发 展 , 也 是 现今 计算 机 技术 发 展 的 结果 . 本 文 讨论 的 是 计算 流体 力学 应 用 问题 , 包括 数学 模型 的 
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选择 , 计算 方法 的 选用 , 及 其 在 飞行 器 设计 中 特殊 问题 上 的 应 用 . 


.学 位 论文 李 永 宏 多 孔 发 汗 冷 却 控制 及 数值 模拟 2007 
多 孔 发 汗 冷却 是 降低 高 温 条 件 下 工作 部 件 (如 火箭 、 导 弹 、 飞 船 等 高 速 飞行 器 的 表面 ) 的 温度 的 有 效 措施 。 为 了 使 高 温 下 工件 材料 表面 不 被 烧 蚀 ， 有 较 长 
的 使 用 寿命 ， 而 又 不 能 改变 材料 的 形状 ， 采 用 多 孔 发 汗 冷 却 技术 是 非常 有 效 的 手段 。 因 此 ， 随 着 工业 技术 的 高 速 发 展 ， 多 孔 发 汗 冷却 技术 已 受到 了 越 来 越 多 
















































































由 于 多 孔 材料 结构 十 分 复杂 ， 所 以 其 内 部 流动 和 换 热 过 程 也 非常 复杂 ， 在 研究 过 程 中 既 要 考虑 消 流 ， 又 要 考虑 传 热 、 传 质 甚 至 化 学 反应 。 内 部 结构 不 规 
律 性 ， 造 成 了 内 部 阻力 的 不 均匀 性 ， 所 以 多 孔 材 料 内 部 的 流动 与 换 热 规律 的 研究 非常 困难 。 目 前 ， 多 孔 发 汗 冷 却 的 大 多 数 研 究 工作 仅 限 于 冷却 剂 对 高 温 气 体 
界 层 的 破坏 以 及 使 对 流 减 弱 的 效果 ， 而 对 多 孔 介质 内 部 的 传 热 过 程 及 温度 分 布 进行 的 研究 很 少 。 对 发 汗 冷却 问题 的 完整 的 物理 模型 和 数学 模型 的 描述 应 该 
括 多 孔 介质 内 部 和 主流 区 的 流动 与 换 热 。 

因此 ， 本 文 将 主要 工作 放 在 多 孔 介 质 内 部 的 传 热 过 程 规律 的 研究 上 。 多 和 孔 介质 由 于 其 内 部 结构 的 复杂 性 和 无 规律 性 ， 使 得 传 热 规律 的 研究 成 为 当前 的 
i 对 复杂 无 规律 的 模型 作 适当 简化 是 首要 的 任务 。 本 文 首先 假定 多 孔 介质 各 向 均匀 同性 ， 在 此 基础 上 推导 出 无 化 学 反应 、 不 可 压缩 发 汗 剂 情况 下 多 和 孔 
部 流体 与 介质 换 热 的 能 量 控 制 方程 ， 然 后 通过 计算 流体 力学 软件 对 多 孔 材 料 发 汗 冷 却 进行 二 维 数值 模拟 计算 。 模 拟 计算 了 不 同 孔 际 率 和 不 同 冷却 剂 流 
况 下 的 多 孔 材 料 内 部 温度 分 布 。 在 温度 分 布 和 速度 分 布 分 析 计 算 的 基础 上 ， 进 一 步 揭 示 了 多 孔 材料 发 汗 冷却 机 理 ， 为 多 孔 发 汗 冷 却 在 工程 应 用 中 的 推广 
































= 















































A 



































































































































7. 会 议论 文 陈 兰 . 胡 汉 东 高 速 旋转 再 入 飞行 器 防 热 层 热 滞 后 数值 计算 2000 
该 文 利 用 数值 方法 研究 了 高 速 飞行 器 有 攻 角 旋 再 入 时 出 现 防 热 层 热 滞后 现象 。 该 方法 在 室 间 离散 上 采用 分 块 的 有 限 体积 法 ， 时 间 推进 采 
塔 积分 。 实 践 表明 该 数值 计算 稳定 ， 时 间 精 度 高 ， 可 以 很 好 的 模拟 复杂 防 热 结构 瞬 态 温度 场 非 定常 热 响应 过 程 。 


8. 会 议论 文 唐 伟 . 桂 业 伟 热 防护 与 热 结 构 性 能 计算 软件 系统 HEATips 简 介 2009 
热 防护 问题 是 制约 高 速 飞行 器 研制 水 平 的 主要 因素 之 一 。 从 计算 网 格 、 热 环境 、 热 响应 、 热 应 力 、 烧 蚀 等 角度 研究 开发 了 热 防护 与 热 结构 性 能 计算 软件 
系统 HEATips， 介 绍 软件 系统 HEATips 的 计算 功能 、 结 构 及 部 分 典型 计算 结果 。 


会 议论 文 WMA. 张 涵 信 . 沈 消 . HE 丙 速 飞行 器 流 场 与 回 体 结 者 构 温 度 场 一 体 化 计算 新 方法 研究 及 其 应 用 2003 

基于 LevelSet 的 原理 , 本 文 探索 和 研究 了 一 种 不 需 作 流 / 固 移 合 迁 代 的 一 体 化 计算 方法 , 通过 将 气体 流动 和 固体 结构 传 热 用 统一 的 方程 组 进行 描述 并 用 统 
-的 方法 进行 求解, 避免 了 流 、 轩 之 间 的 耦合 迭代 问题 .为 了 验证 方法 的 可 靠 性 , 本 文 对 高 烩 高 超声 速 气流 绕 二 维 不 锈 钢 圆 管 的 流动 、 气 动 加 热 以 及 结构 热 响 
应 问题 进行 了 计算 , 结果 与 实验 符合 很 好 . 在 此 基础 上 本 文 将 方法 成 功 应 用 到 了 高 速 飞机 沿 飞行 轨道 飞行 时 的 热 环境 和 结构 热 响应 计算 中 


0. 学 位 论文 EA 高 速 清流 流 场 气动 光学 传输 效应 的 计算 机 仿真 研究 2007 
随 着 新 军事 技术 变革 的 到 来 ， 各 种 无 人 驾驶 精确 制导 的 飞行 器 对 信息 获取 的 精确 性 和 时 效 性 的 依赖 程度 大 幅度 地 提高 ， 采 用 光学 成 像 探 测 跟踪 制导 和 追 
求 高 速 飞 行 ， 已 成 为 当代 各 种 精确 制导 武器 的 必然 发 展 趋势 。 随 着 这 两 项 技术 的 广泛 应 用 ， 气 动 光 学 问题 也 应 运 而 生 。 飞 行 器 速度 越 高 ， 采 用 光学 成 像 探测 
制导 系统 的 精确 制导 武器 气动 光学 效应 越 尽 重 ， 并 严重 影响 到 探测 的 精度 。 气 动 光 学 问题 已 经 成 为 这 两 种 技术 同时 应 用 的 瓶 项 ， 该 领域 的 研究 进展 将 直接 影 
响 高 速 飞行 器 光学 探测 系统 的 精确 性 ， 影 响 新 一 代 精 确 打击 武器 的 发 展 。 气 动 光 学 效应 包括 : 气动 热效应 、 气 动 热 辐 射 效应 和 气动 光学 传输 效应 。 前 两 种 效 
降低 了 导 引 头 对 目标 的 探测 的 信 品 比 ， 后 一 种 效应 降低 了 导 引 头 对 目标 的 探测 、 跟 踪 和 识别 的 能 
气动 光学 所 处 的 工作 条 件 相当 苛刻 ， 除 了 超 音速 气流 的 要 求 外 ， 往 往 还 涉及 高 温 、 高 压 、 强 光 辐 射 、 有 毒 和 强 腐蚀 等 复杂 环境 ， 试 验 成 本 相当 高 ， 因 而 
数值 模拟 已 经 成 为 气动 光学 理论 研究 的 一 个 重要 手段 和 工具 。 本 研究 就 是 针对 气动 光学 传输 效应 的 基础 性 问题 : 流 场 介质 折射 率 变化 是 影响 光学 传输 最 直接 

的 因素 ， 即 针对 光线 穿 过 高 速 消 流 流 场 作 定性 分 析 和 定量 计算 而 开展 的 气动 光学 传输 效应 计算 机 仿真 研究 工作 。 
本 文 在 概述 了 气动 光学 效应 的 全 貌 基础 上 ， 研 究 了 国内 外 气动 光学 效应 及 计算 机 仿真 的 历史 与 现状 。 在 描述 气动 光学 传输 效应 光学 理论 之 后 ， 分 析 了 气 
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动 光学 传输 效应 定量 研究 方法 ， 


结合 现 有 的 高 速 清流 流 场 数据 ， 
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了 基于 随机 相位 屏 的 相位 屏 思想 与 基于 消 流 结构 研究 的 分 界面 厚度 方法 ， 研 





究 了 分 界面 























厚度 方法 在 气动 光学 传输 效应 计算 机 仿真 中 的 作用 ， 并 实现 了 一 个 高 速 流 场 气动 光学 传输 效应 的 计算 机 仿真 软件 。 

在 课题 小 组 的 研究 的 基础 上 ， 着 重 进行 了 如 下 的 研究 

Le 气动 光学 传输 效应 的 定量 研究 方法 气动 光学 的 定量 研究 是 成 功 运用 计算 机 进行 仿真 的 关键 ， 方 法 的 选择 对 于 仿真 结果 有 着 重要 的 影响 。 因 此 本 文 从 
工程 应 用 角度 研究 了 气动 光学 传输 效应 的 定量 研究 方法 ， 即 传统 的 实验 方法 、 统 计 方法 和 近 些 年 兴起 的 随机 相位 屏 法 、 满 流 结构 研究 方法 。 并 利用 随机 相位 














屏 方法 实现 了 大 气 湛 流 的 仿 

2. 分 析 了 现 有 
文 分 别 从 概率 、 统 计 、 
分 界面 厚度 方法 进行 了 气动 
位 屏 法 和 折射 率 分 界面 方法 针对 4 


3. 利 


结合 


本 文中 


本 文 链 





随机 相 
场 气动 光学 传输 效应 计算 机 仿真 软 
究 表明 : 高 速 清流 
效应 进行 仿真 ， 仿 真光 程 差 的 相对 误差 4 
动 光学 失真 








分 布 函数 、 











[计算 ， 为 后 续 
高速 湛 流 流 场 数据 














件 的 设计 与 实现 。 





的 仿真 研 








的 研究 。 





的 气动 光学 仿真 打下 了 基础 。 
分 析 现 有 CFD (计算 流体 力学 ) 流 
图 像 重 现 等 多 角度 分 析 了 、 比 较 了 流 场 数据 。 
光学 失真 波 前 仿真 详细 介绍 了 基 了 
特定 速度 ， 高 度 的 流 场 进行 了 气动 光学 失真 波 前 





Bai» ATRE 





FP 高 梯度 分 界面 产生 大 尺度 气动 光学 失真 ， 能 较 好 利 月 
下 超过 4% 。 分 界面 厚度 方法 对 于 气动 光学 失真 的 预测 











E 解 流 场 特征 ， 合 理 正确 采 


仿真 方法 有 着 








FF 满 流 结构 研究 的 分 界面 厚度 模型 ， 并 通过 ei 


的 仿真 ， 分 析 了 结果 。 4. 

















接 : http://d. g. wanfangdata. com. cn/Thesis Y1102098. aspx 
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onal 方程 推导 了 分 界面 
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- 定 的 指导 意义 。 本 





厚度 方法 。 
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折射 率 分 界面 厚度 模拟 气动 光学 效应 ， 采 
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